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ABSTRACT 


Recent  studies  of  the  mechanism  involved  when  a  reservoir  fluid 
is  displaced  miscibly  by  a  light  hydrocarbon  mixture  have  resulted  in  the 
development  of  theoretical  methods  to  predict  the  approximate  conditions 
necessary  in  order  that  miscibility  may  occur.  One  of  these  methods  employs 
the  use  of  a  pseudo -ternary  phase  diagram  and  requires  the  calculation  of 
the  critical  temperature  and  pressure  of  a  mixture  of  both  the  displacing 
and  displaced  phases.  A  brief  theoretical  analysis  of  the  procedure  was 
undertaken  and  comparisons  were  made  between  the  method  used  to  calculate 
criticals  by  the  prediction  correlation  and  newer  methods  of  critical  calcu¬ 
lation. 

The  second  part  of  the  study  is  devoted  to  examining  the  mechanism 
of  miscible  fluid  displacement  in  liquid  systems  and  how  it  is  influenced 
by  viscous  fingering.  A  unique  method  of  recording  pressure  profiles  as  a 
function  of  time  during  an  actual  miscible  displacement  has  been  devised. 
Production  data  and  pressure  profiles  along  the  core  for  various  time 
increments  were  obtained  and  are  plotted  and  discussed.  Experimental  results 
indicate  that  ideal  miscible  displacements  take  place  only  at  low  viscosity 
ratios  and  that  laboratory  verification  of  theoretical  miscible  predictions 
requires  strict  control  of  all  causes  of  poor  core  sweep  efficiency. 
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UTRODUCTION 

The  ideal  miscible  displacement  of  reservoir  oil  depends  upon 
the  formation  of  a  graded  staturation  bank  between  the  original  oil  in  place 
and  the  injected  displacing  fluid.  This  bank  is  all  in  a  single  phase  with 
no  interfaces  present  and  therefore  there  are  no  capillary  forces  to  cause  the 
oil  to  be  retained  in  pore  channels  as  residual  oil.  This  situation  can  be 
achieved  by  employing  one  or  more  different  methods^  namely  LPG  slug  process, 
solvent  flooding,  high  pressure  gas  drive  or  condensing  gas  drive.  These 
processes  are  briefly  outlined  below?- 

1)  Miscible  LPG  Slug:  A  predetermined  quantity  of  liquified  petroleum 
gas,  usually  propane,  is  injected  into  the  reservoir  initially  and  is  driven 
toward  the  producing  wells  by  a  secondary  phase  of  water  or  dry  gas. 

2)  Solvent  Flooding:  The  entire  reservoir  is  flooded  with  a  solvent 
such  as  alcohol  or  LPG  and  some  provision  is  made  for  recovering  the  solvent 
at  the  end  of  the  flood. 

3)  High  Pressure  Gas  Drive:  Lean  gas  is  injected  into  the  reservoir  at 
very  high  pressures  and  this  injected  gas  absorbs  intermediate  hydrocarbons 
(c2  "  Cg)  from  the  reservoir  oil.  This  enrichment  continues  until  a  graded 
saturation  bank  is  built  up  at  the  flood  front  thereafter  miscibly  displacing 
the  reservoir  oil  ahead  of  it. 

4)  Enriched  or  Condensing  Gas  Drive:  In  this  process  gas  rich  in 
Intermediate  hydrocarbons  (Cg  -  C^)  is  injected  into  the  reservoir  under 
conditions  that  will  allow  a  large  portion  of  these  intermediate  components 

to  condense  upon  contact  with  the  reservoir  oil.  Thus  a  low  viscosity,  graded 
saturation  bank  is  formed  which  is  driven  toward  the  producing  wells  by  the 
subsequently  Injected  uncondensed  gas. 
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Three  of  the  factors  which  influence  the  efficiency  of  miscible 
drives  are  the  reservoir  temperature  and  pressure  and  the  composition  of  the 
injected  fluid.  A  certain  pressure  for  a  given  temperature  must  prevail  to 
maintain  miscible  conditions  between  oil  and  injected  fluid.  Proper  require¬ 
ments  for  pressure  control  of  the  reservoir  must  therefore  be  known  in  advance 
of  commencement  of  a  project  for  properly  controlled  miscible  drive  operations. 
The  intermediate  composition  of  the  injected  gas  must  also  be  controlled  such 
that  the  composition  will  be  a  critical  one  and  miscible  drive  can  be 
achieved . 

The  mechanisms  of  high  pressure  gas  drive  and  condensing  gas 

drive  have  been  described  for  the  simple  case  of  three  pure  components. 

(3  33) 

Certain  authors  3  have  endeavoured  to  use  this  description  to  determine 
quantitatively  the  conditions  under  which  miscibility  will  be  obtained  for 
the  case  of  crude  oil  and  natural  gas.  This  requires  that  the  several  com¬ 
ponents  of  the  fluids  be  divided  into  three  groups ,  each  of  which  behaves  like 
a  pure  substance. 

The  primary  objectives  of  this  research  project  were  to  carry  out 
preliminary  theoretical  investigations  of  prediction  of  miscibility  and  to 
experimentally  study  the  mechanism  of  a  miscible  displacement.  It  should 
be  pointed  out  that  the  test  results  reported  are  for  a  liquid-liquid  miscible 
drive 9  where  it  is  known,  miscibility  exists.  This  system  was  selected  so 
that  quantitative  relationships  between  viscous  fingering  and  viscosity  ratio 


could  be  obtained. 
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PREDICTION  OF  MISCIBILITY 

Theory  of  Benham,  Dowden  and  Kunzman. 

A  recent  analysis  of  the  mechanism  of  miscible  displacement  by 
Benham,  Dowden  and  Kunzman,  has  led  to  the  development  of  a  simplified 
method  for  predicting  conditions  of  miscibility  for  an  enriched  gas  drive. 

Their  work  makes  use  of  a  pseudo-ternary  diagram.  Figure  1,  in  which  the 
mixture  of  reservoir  oil  and  enriched  injected  gas  is  divided  into  three 
components,  namely  and  C£  -  from  the  injected  gas  and  C^+  from  the 
reservoir  oil.  A  vertical  limiting  tie  line  or  limiting  tie  line  of 
infinite  slope  is  then  imposed  upon  the  diagram  to  indicate  the  phase 
conditions  which  can  prevail  during  miscible  displacement. 

Justification  of  the  use  of  this  type  of  ternary  representation 
is  reported  as  follows:  "if  one  considers  a  position  some  distance  from  the 
injection  point,  the  character  of  the  C^+  fraction  in  the  liquid  phase  probably 
will  not  change  by  an  appreciable  amount  during  its  approach  to  miscibility. 
Also  if  considerations  are  limited  to  the  region  close  to  the  critical  point 
and  the  associated  limiting  tie  line,  the  character  of  the  intermediate 
fraction  will  be  close  to  that  of  the  injection  mixture.  Therefore,  a  pseudo¬ 
ternary  phase  diagram  (drawn  using  intermediates  of  character  equal  to  that 
of  the  injection  mixture  and  C^+  fractions  having  a  character  equal  to  that 
of  the  reservoir  fluid)  will  have  quantitative  significance  in  the  region  of 
the  critical  point". 

The  use  of  the  vertical  limiting  tie-line  concept  is  explained 
by  the  fact  that  most  hydrocarbon  systems  which  were  studied  exhibited 
limiting  tie  lines  of  negative  slope.  Applying  the  afore  stated  assumption 
it  is  then  concluded  that:  "the  limiting,  miscible,  injection  fluid 


-  4  - 


FIG.  NO.  I 

MISCIBLE  DISPLACEMENT  OF  RESERVOIR  FLUID 

BY 

RICH  GAS 

(AFTER  BENHAM,  DOWDEN  AND  KUNZMAN) 


MOL  % 


composition  may  be  calculated  by  combining  reservoir  fluid  with  various 
mixtures  of  potential  displacing  fluid  LFG  and  methane  in  various  proportions 
and  calculating  the  critical  temperature  and  pressure  for  those  mixtures". 


A  number  of  such  calculations  were  carried  out  and  the  results 
for  one  particular  reservoir  fluid  and  a  number  of  different  displacing 
fluid  mixtures  are  shown  on  Figure  2.  Ey  extending  these  calculations  to 
cover  a  number  of  different  reservoir  fluids  and  displacing  fluid  compositions, 
a  generalized  correlation  was  developed.  The  final  result  of  this  work  was 
a  correlation  of  the  maximum  permissible  methane  concentration  in  the  displacing 
fluid  for  a  miscible  displacement  of  reservoir  fluids  as  a  function  of 
reservoir  temperature,  reservoir  pressure,  CL  molecular  weight  of  the  reservoir 
fluid  and  C2+  molecular  weight  of  the  displacing  fluid.  An  example  of  one 
of  the  final  charts  is  shown  on  Figure  3 • 

The  use  of  a  chart  such  as  this  entails 

1)  Calculation  of  the  molecular  weight  of  the  C^+  fraction  of  the 
reservoir  fluid. 

2)  Calculation  of  the  molecular  weight  of  the  C2  -  C4  fraction  of 
the  displacing  fluid. 

3)  Entry  onto  the  correct  correlation  chart  for  the  reservoir 
pressure  and  temperature  in  question  and  subsequent  determination 
of  the  maximum  permissible  methane  concentration  in  the 


displacing  fluid. 


C 


" 
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FIG.  NO.  2 

CRITICAL  POINTS  FOR  MIXTURES 

OF 

RESERVOIR  FLUID  E,  PROPANE  AND  METHANE 
(AFTER  BENHAM,  DOWDEN  AND  KUNZMAN) 
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Discussion  of  Benham,  Dowden  and  Kunzman  Correlation. 

Work  recently  undertaken  at  the  University  of  Alberta  has  led 
to  some  interesting  speculation  as  to  the  feasibility  of  predicting  miscible 
conditions,  whether  it  be  by  the  method  of  Benham,  Dowden  and  Kunzman 
or  by  some  other  procedure.  Use  of  the  pseudo -ternary  type  representation 
for  a  multicomponent  mixture  has  been  questioned  by  recent  authors^)  on 
the  basis  that  it  does  not  represent  the  true  composition  of  various  mixtures 
of  gas  and  liquid  in  each  section  of  the  porous  media.  It  is  contended  that 
the  composition  of  any  mixture  in  a  particular  section  of  the  porous  media 
depends  upon  the  nature  of  the  porous  media  and  the  velocity  of  displacement. 

The  mechanism  by  which  a  wet  gas  displaces  a  crude  oil  (condensing 
gas  drive)  is,  for  a  period  of  time,  a  non-equilibrium  gas  displacement 
process.  In  the  actual  displacement  process  a  saturation  profile  (gas-oil) 
in  terms  of  distance  is  initially  established  and  in  some  cases  may  continue 
to  grow. 

Even  supposing  that  the  pseudo -ternary  representation  can  be 
considered  thermodynamically  rigorous  in  a  porous  media,  certain  other  aspects 
of  it  must  be  considered.  Benham  et  aX  assumed  that  the  limiting  tie  line 
would  be  of  negative  slope  and  therefore  that  a  vertical  limiting  tie  line 
would  give  a  conservative  estimate  as  to  the  amount  of  methane  to  be  used 
for  injection  fluid  dilution.  This  assumption  can  probably  be  considered 
correct  provided  the  amount  of  ethane  in  the  injection  fluid  is  not  excessive. 
No  consideration,  however,  was  taken  for  the  presence  of  non-hydrocarbon 
components  (H2S,  C02,  N2)  in  the  mixture.  The  presence  of  such  components 
will  undoubtedly  affect  the  critical  temperature  and  pressure  of  the  mixture 
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as  well  as  the  shape  of  the  two-phase  region  on  the  diagram,.  Whether  or  not 
the  assumption  of  a  vertical  limiting  tie  line  will  hold  in  cases  of  mixtures 
containing  appreciable  amounts  of  non-hydrocarbon  components  should  be 
determined.  Recent  studies ^33)  have  indicated  that  different  groupings  of 
components  can  be  used  when  constructing  a  ternary  diagram.  It  has  been 
suggested  that  the  best  ternary  representation  can  be  obtained  by  dividing 
a  mixture  into  the  components 

1)  C^,  CO2 ,  and  N2 

2)  C2-C6  and  H2S 

3)  <^+ 

Critical  conditions  for  mixtures  studied  by  Benham,  Dowden  and 

( ig ) 

Kunzman  were  calculated  using  the  method  of  Kurata  and  Katzv  y ’  as  modified 
by  Davis,  Bertuzzi,  Gore  and  Kurata^.  Various  mixtures  of  displacing 
fluid  and  reservoir  fluid  were  used  to  calculate  these  criticals.  It  is 
therefore  evident  that  the  correct  proportions  of  reservoir  fluid  and 
displacing  fluid  along  with  an  accurate  method  of  calculating  criticals 
must  be  employed  in  order  that  the  final  correlations  have  some  degree  of 
accuracy. 

Knowing  the  temperature  and  pressure  of  the  reservoir  and  knowing 
the  composition  of  the  reservoir  fluid,  one  can,  by  a  series  of  trial  and 
error  calculations,  determine  a  mixture  which  will  have  as  its  critical 
properties  the  pressure  and  temperature  conditions  of  the  reservoir.  Once 
the  composition  of  this  mixture  has  been  obtained  it  is  possible  to  ascertain 
the  proportions  of  displacing  fluid  and  reservoir  fluid  which  are  in 
equilibrium  at  this  condition.  However,  can  it  be  assumed  that  the 
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proportions  of  these  fluids  will  remain  the  same  throughout  the  critical 
zone  in  the  porous  media?  If  they  do  not  then  loss  of  miscibility  can 
occur. 
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Alternate  Theoretical  Considerations 

It  was  originally  thought  that  the  controlling  factor  in 

determining  the  accuracy  of  the  Benham  et  al.  prediction  was  the  calculation 

of  the  critical  conditions.  With  this  in  mind  an  investigation  of  alternate 

methods  of  calculating  criticals  was  undertaken.  It  was  found  that  many 

methods  have  been  developed  to  predict  the  critical  temperature  and  pressure 

of  mixtures  and  they  are  presented  in  literature  10*12*19*22*25*29)  ^ 

(29) 

earliest  of  these  was  in  1937  when  Smith  and  Watson"  developed  a 

correlation  to  be  used  on  petroleum  fractions.  Their  method  requires  a 

knowledge  of  the  gravity  of  the  stock*  of  its  molecular  weight*  and  of 

the  molal  average  and  weighted  average  boiling  points. 

(22) 

Mayfield  ’  in  1942  proposed  a  method  for  binary  paraffinic 

systems  which  makes  use  of  a  weighted  average  critical  temperature.  A 

critical  pressure  is  calculated  using  the  differences  between  the  weighted 

average  critical  pressure  and  the  experimental  literature  values. 

In  1942*  Kurata  and  Katz ^9)  presented  a  method  which  is  an 

extension  of  Smith  and  Watson's  technique.  Their  calculation  uses  molecular 

weight*  pseudo-critical  temperature  and  a  new  function*  M*  designed  to 

represent  boiling  range. 

(25) 

Organ! ck  •  1  in  1953  was  the  first  to  develop  a  calculation 

specifically  for  a  multi -component  system  containing  hydrogen.  This  method 
uses  a  weighted  average  for  calculation  of  critical  temperature.  To  calculate 
critical  pressure  by  this  method  however  it  is  necessary  to  know  the  molal 
average  boiling  point  and  the  weight  average  equivalent  molecular  weight. 

In.  1953 )  Davis*  Bertuzzl*  Gore  and  Kurata ^ ®  ^  modified  the 
method  which  bad  previously  been  developed  by  Kurata  and  Katz.  They 
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corrected  the  weighted  average  critical  temperature  for  each  binary  pair 
in  the  system.  This  was  done  because  deviation  is  greatest  at  nearly  equal 
weight  fractions.  They  also  modified  and  substituted  an  arbitrary  variable, 
M*,  for  the  M  function  in  the  Kurata  and  Katz  method.  This  was  done  to 
correct  for  the  presence  of  nitrogen. 

Eilerts  in  1957  developed  a  series  of  correlations  for 
binary  and  complex  mixtures  based  upon  the  boiling  points  of  the  pure 
components,  the  molal  average  boiling  point,  the  pseudo-critical  pressure, 
and  a  weighted  combination  of  the  critical  temperatures  and  boiling  points 
of  the  pure  components. 

In  i960.  Grieves  and  Thodos'-^)  presented  a  new  approach  for 
calculating  critical  temperatures  and  pressures  of  mixtures  containing 
aliphatic,  napathenic  and  aromatic  hydrocarbons  as  well  as  hydrogen.  They 
introduced  two  dimensionless  temperature  parameters  (V  and  9)  which  are 
defined  by  the  molal  average  boiling  point,  the  boiling  point  and  the  dew 
point  of  the  mixture,  all  at  atmospheric  pressure.  This  method  was 
developed  primarily  for  binary  systems  and  consequently  it  has  not  been  used 
to  any  great  extent  on  complex  mixtures. 

After  a  preliminary  examination  of  the  various  methods  of 
calculating  criticals  was  completed  it  was  decided  to  examine  three  of  them 
in  detail.  The  methods  which  were  deemed  to  be  best  suitable  for  calculating 
critical  conditions  of  a  miscible  type  mixture  were  Organick's;  the  Davis, 
Bertuzzi,  Gore  and  Kurata  modification  and  Eilerts'. 
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Method  of  Organick. 

Critical  Temperature  and  Pressure . 

Organick* s  method  employs  the  use  of  two  properties 

1)  Molal  average  boiling  point. 

2)  Weight  average  equivalent  molecular  weight. 

The  molal  average  boiling  point  is  calculated  by  the  relationship 

B  =  n^ 

where  B  =  molal  average  boiling  point 

=  mole  fraction  of  component  "i" 

Bj_  =  normal  boiling  point  for  component  "i". 

To  obtain  the  molal  average  boiling  point  of  the  Cy+  fraction, 
Organick  used  a  correlation  based  on  the  A.S.T.M.  distillation  curve  of  the 
C y+  fraction.  This  correlation  relates  the  slope  of  the  10$  to  90$  portion 
of  the  curve  to  the  volumetric  boiling  point  and  the  true  boiling  point. 

Using  this  correlation  it  is  possible  to  obtain  the  mole  fraction  as  well  as 
the  molal  average  boiling  point  of  the  Cy+  fraction. 

The  "weight  average  equivalent  molecular  weight"  is  determined 
by  the  relationship 


W  =  mi  Mjgi 

where  W  =  weight  average  equivalent,  molecular  weight 
m^  =  weight  fraction  of  component  "i" 

MEi  =  equivalent  molecular  weight  of  component  "i 
i)  for  paraffin  hydrocarbons.  Mg  =  M 
ii)  for  Isomers  or  olefins,  Mg  is  the  molecular  weight 

that  a  paraffin  would  have  if  it  had  the  same  boiling 
point  as  the  isomer  or  olefin. 


. 


iii)  for  aromatic  or  napthenic  compounds  a  special 

characterization  factor  is  used. 

Organick  then  plotted  Pc  vs  T0  for  various  mixtures  and  using 

constant  B's  as  a  parameter,  drew  in  families  of  lines.  In  a  like  manner 

he  plotted  P  vs  T  again  and  using  constant  W?s  as  parameters,  drew  in 
c  c 

another  set  of  lines.  He  then  superimposed  one  plot  upon  the  other  and 
obtained  a  plot  of  Pc  vs  Tc  with  a  grid  of  constant  B’s  and  W’s  upon  it. 
Organick1 s  investigations  were  carried  out  in  two  specific  areas t- 

1)  Mixtures  with  B8s  from  700°R  to  900°R  and  VPs 
from  120  to  220. 

2)  Mixtures  with  B* s  from  250°R  to  360°R  and  W’s 
from  35  to  120. 

In  the  final  correlation  Organick  joined  these  two  plots 
together  by  interpolating  between  known  values.  Therefore  on  his  final 
chart  the  whole  center  section  has  been  dotted  in. 

The  use  of  the  Organick  method  simply  entails  the  calculation 
of  B  and  W  for  a  mixture  and  then  entering  the  chart  to  pick  off  the 
corresponding  Pc  and  Tc  for  that  particular  B  and  W. 

Method  of  Kurata  and  Katz ;  -  as  modified  by  Davis  et.al. 

Critical  Temperatures. 

The  critical  temperatures  of  many  mixtures  are  found  to  be 
nearly  equal  to  the  weight  average  critical  temperatures,  mTc, ,  of  these 
mixtures.  Kurata  and  Katz  correlated  the  critical  temperatures  of 
hydrocarbon  mixtures  using  the  molal  average  critical  temperature  ’  Tcn> 
and  the  average  molecular  weight,  with  a  special  treatment  for  gases 


containing  methane . 
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The  method  developed  by  Davis  et  al.  for  volatile  mixtures 
consists  of  applying  a  set  of  corrections  to  mTc,  the  weight  average 
critical  temperature.  Existing  data  shows  that  mTc  approximates  the  true 
critical  temperature,  T  .  The  deviation  is  greatest  at  nearly-equal  weight 
fractions  of  any  two  components  and  can  be  approximated  by : - 

Tc  -  mTc  =  A^gmpmg 

where  m.  and  m2  are  the  weight  fractions  of  components  1  and  2  respectively. 
The  constant  A-^2  is  characteristic  of  the  particular  binary  mixture.  This 
relationship  was  extended  to  complex  mixtures  by  the  equation:- 

Tc  -  mTc  “  Al,2mlm2  Al,  3mlm3 . 

A2,3m2®3  . 

where  a  term  is  included  for  every  possible  binary  pair. 

Davis  et  al.  derived  constants  for  the  above  equation  from  existing 
data  on  binary  systems.  Where  existing  data  were  not  available,  interpolation 
and  extrapolation  procedures  were  used.  A  number  of  constants  were  also 
obtained  from  data  on  complex  mixtures  but,  according  to  the  authors,  are 
subject  to  considerable  error. 

Critical  Pressure 

The  critical  pressures  of  the  mixtures  studied  by  Davis  et  al.  were 

correlated  by  a  modification  of  the  method  of  Kurata  and  Katz.  The  original 

correlation  was  a  plot  of  TCP  cq/Pc  Vs  T  cn  with  M  as  a  parameter.  This 

o 

correlation  was  modified  in  the  low  temperature  region  (T  cn  400  R)  by 
substituting  an  arbitrary  variable,  M*,  for  molecular  weight,  M,  in  order 
to  account  for  gases  containing  nitrogen.  When  no  nitrogen  is  present,  M* 


is  equal  to  M. 
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Method  of  Eilerts s- 

Critlcal  Temperature, 

Eilerts J  method  is  based  upon  the  fact  that  the  critical 
temperature  of  a  mixture  is  higher  than  the  critical  temperature  of  the 
lightest  component  and  lower  than  the  critical  temperature  of  the  heaviest 
component o  Therefore ,  by  averaging  the  critical  temperatures  of  the 
components  of  a  multicomponent  mixture  on  the  basis  of  moial  concentrate  on 
and  giving  additional  weight  to  the  heaviest  components,  he  proposed  that  a 
critical  temperature  can  he  estimated  that  will  approximate  the  measured 
critical  temperature.  The  normal  absolute  boiling  point  of  the  respective 
components  was  used  as  a  weighting  factor  in  order  to  emphasize  the  effect 
of  the  critical  temperature  of  the  components  with  high  molecular  weights 
on  the  estimated  critical  temperature. 

Using  this  principle,  Eilerts  then  obtained  a  correlation  to 
improve  the  accuracy  in  computing  the  critical  temperature  of  hydrocarbon 
mixtures.  The  ratio  of  the  measured  critical  temperature  to  estimated 
critical  temperature ,  Tern /'Ice ,  was  used  in  preparing  a  correlation  in  which 
the  measured  critical  temperature  was  related  to  the  function,  (%  -  Bj/i/Bp 
where:  "  %  =  normal  boiling  point  of  heavy  (liquid) 

portion  of  mixture,  °R 

Bp  =  normal  boiling  point  of  light  (vapor)  portion 
of  mixture,  °R 

Bp  =  normal  boiling  point  of  entire  mixture,  °R. 

The  correlating  function  depends  for  value  solely  on  the  composition  of  the 
mixture.  The  final  Eilerts’  corre.lat.ion  plotted  the  ratio  of  the  measured 
critical  temperature  to  estimated  critical  temperature  vs  a  composition 


. 
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variable,  Ns,  with  selected  values  of  Bg/Bg  as  parameters.,  The 
composition  variable,  !NS,  was  designated  as  the  mole  fraction  of  the 
mixture  which  is  in  the  gas  phase  at  l4„7  psia  and  75°F,  The  values  for 
Tcm  and  Ns  were  obtained  from  experimental  data  while  TC0,  Bg  and  Bg  were 
calculated . 

Critical  Pressure . 

Eilerts  realized  that  the  critical  pressure  of  a  natural 
hydrocarbon  mixture  is  usually  higher  than  that  of  any  component  in  the 
mixture  and  also  higher  than  the  molal  average  of  the  critical  pressures  of 
all  components.  Although  computation  of  the  molal-average  critical  pressure 
can  be  an  important  step  in  computing  the  true  critical  pressure  of  a  mixture, 
the  molal-average  critical  pressure  alone  is  not  a  reliable  indication  of  the 
true  critical  pressure.  Eilerts  developed  a  correlation  plotting  the  ratio 
of  the  measured  critical  pressure  to  the  molal-average  critical  pressure  vs 
Ks  with  Bg/ Bg  as  parameters.  Again  Pcm,  the  measured  critical  pressure,  and 
Ws,  the  mole  fraction  of  the  mixture  in  the  vapor  phase  at  14.7  psia  and  75°F, 
were  obtained  from  experimental  data.  The  molal-average  critical  pressure, 
Pen*  and  the  boiling  points  of  the  liquid  and  vapor  portions  of  the  mixture, 

Bg  and  Bg  respectively,  were  calculated. 

Use  of  the  Eilerts*  correlations  for  determining  critical 
temperatures  and  pressures  of  mixtures  entails  calculation  of  the  vapor- 
liquid  equilibrium  ratios  of  each  component  which  would  occur  at  3,4.7  psia 
and  75°F.  This,  of  course,  requires  an  accurate  analysis  of  the  components 
in  the  Cy +  portion  of  the  mixture.  To  obtain  this  analysis,  Eilerts  uses 
the  results  from  a  U.S.B.M.  distillation  on  the  Cy+  portion  of  the  mixhu 
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A  comparison  of  the  critical  temperatures  and  pressures  obtained 
using  the  three  different  calculation  methods  on  various  naturally  occurring 
hydrocarbon  mixtures  is  presented  in  Table  I  of  the  Discussion  and  Results 
Section. 
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Experimental  Apparatus: 

Figure  4  is  a  schematic  diagram  of  the  experimental  apparatus  used 
during  this  project.  The  simulated  reservoir  employed  was  a  consolidated  Berea 
sandstone  core,  2  inches  in  diameter  and  3 6  inches  long.  The  core  was 
mounted  in  a  length  of  steel  casing  4.5  inches  in  diameter.  The  annular 
space  between  the  casing  and  the  core  was  filled  with  Shell  Epon  Resin.  A 
continuous  bond  between  the  core  and  the  resin  and  between  the  resin  and  the 
casing  wall  was  thereby  obtained.  This  insured  that  there  would  be  no  bypassing 
of  fluids  past  the  core.  The  average  porosity  of  the  core  was  experimentally 
determined  to  be  16.3  per  cent  while  the  average  liquid  permeability  was  266 
millidarcies.  The  casing  and  core  were  mounted  between  two  flanges  which  were 
held  together  with  tie -rods. 

Pressure  taps  were  drilled  at  6  inch  intervals  along  the  core  and 

Statham  Model  PM80Tct  15-350  Differential  Pressure  Transducers  were  mounted 

in  series  along  the  core.  Originally  four  of  these  transducers  were  used  but 

later  this  number  was  changed  to  three  in  order  that  the  sections  of  the  core, 

across  which  each  transducer  measured  the  differential  pressure,  would  be  even 

in  length.  The  transducers  had  a  pressure  range  of  plus  or  minus  15  pounds 

per  square  inch  differential  and  an  absolute  operating  pressure  of  5000  p.s.i.a. 

o  o 

They  were  temperature  compensated  in  the  interval  0  F  to  200  F.  The  maximum 
excitation  voltage  recommended  for  each  transducer  was  16  volts  and  the 
approximate  output  resistance  of  each  was  350  ohms.  Input  voltage  to  the 
transducers  was  supplied  by  a  Kepco  Model  SC  -  18  -  1,  voltage  regulated  D.C. 
power  supply.  The  D.C.  output  of  this  instrument  was  0  to  18  volts  and  0  to 
1  amp.  with  an  A.C.  input  of  105  to  125  volts  at  50  to  65  cps.  All  voltages 
within  the  system  were  measured  with  a  Non-Linear  Systems  Model  484a  Industrial 
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Digital  Voltmeter. 

The  output  from  the  pressure  transducers  was  recorded  'by  a.  Honey- 
well  Model  906c  Visicorder.  This  instrument  is  a  direct  reading  type  of 
oscillograph  which  prints  a  trace  on  continuously  moving  photographic  paper. 

This  is  done  ty  passing  the  input  current  from  each  transducer  to  a  miniature 
mirror  galvanometer  in  the  Visicorder.  The  mirror  galvanometer  then  reflects 
a  beam  of  ultra-violet  light  onto  a  moving  sheet  of  photographic  paper.  Using 
this  equipment  a  permanent  record  of  differential  pressure  versus  time  between 
any  two  pressure  taps  could  he  obtained  throughout  a  complete  displacement. 
Figure  5  is  a-  reproduction  of  a  chart  obtained  from  the  Visicorder  during  an 
actual  displacement. 

The  pressure  sensing  and  recording  system  was  calibrated  prior 
to  the  initial  displacement  of  fluid  into  the  core.  This  was  done  by  placing 
a  full-scale  differential  pressure  across  each  transducer  and  varying  the 
input  voltage  to  the  transducer  until  a  full-scale  reading  was  obtained  on 
the  Visicorder.  The  differential  pressure  on  the  transducer  was  then  lowered 
in  increments  to  zero  and  increased  to  full-scale  to  check  the  Visicorder  for 
linearity.  All  pressure  differentials  were  measured  by  a  mercury  manometer 
during  this  calibration  procedure.  An  attempt  was  made  to  check  the  calibration 
factor  of  each  transducer  as  reported  by  the  manufacturer  hut  this  became 
impossible  due  to  the  insensitivity  of  voltage  measuring  devices  which  were 
available.  It  was,  however,  discovered  that  full-scale  deflection  of  the 
diaphragm  in  the  transducer  required  approximately  0.1  cc  of  liquid. 

The  core  and  transducers  were  enclosed  in  a  constant  temperature 
air  hath  while  the  recording  instruments  were  mounted  in  an  adjacent  location. 
Displacements  were  carried  out  using  a  Ruska  Model  2243  -  Bill  -  U  Q  Positive 
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Displacement  Pump.  This  pump  was  of  the  double -cylinder  type  and  was  geared 
for  discharge  rates  of  from  10  to  1200  cc  per  hour  per  cylinder  and  calibrated 
to  0.01  cc.  The  pump  was  manifolded  to  a  series  of  high  pressure  sample 
containers  such  that  simultaneous  injection  of  two  different  fluids  could  be 
carried  out  if  required.  All  absolute  displacement  pressures  as  well  as  the 
upstream  and  downstream  pressure  on  the  core  were  measured  with  a  3000  p.s.i.a. 
Heise  gauge.  A  Grove  back  pressure  regulator,  supplied  by  a  nitrogen  bottle, 
was  used  to  maintain  a  back  pressure  on  the  core  during  displacements. 

The  effluent  produced  during  a  displacement  was  ran  through  a 
Ruska  Model  7598  Flash  Separator.  This  separator  had  a  pressure  range  of  0 
to  500  p.s.i.g.  and  a  graduated  glass  separation  column  with  a  volume  of  50  cc. 
Any  gas  produced  from  the  separator  was  measured  with  a  wet  test  meter  before 
being  vented  to  the  atmosphere. 
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Experimental  Procedure  ; 

Prior  to  the  introduction  of  any  liquid  into  the  core  the  entire 
apparatus  was  pressure  tested  to  2000  p.s.i.g.  with  nitrogen  to  check  for 
leaks.  A  vacuum  was  also  imposed  upon  the  core  for  2b  hours  as  an  additional 
check.  Air  was  then  flowed  through  the  core  at  various  rates  in  order  to 
obtain  the  permeability  of  the  core  to  air  as  well  as  to  check  the  operation 
of  the  pressure  sensing  and  recording  devices.  The  data  obtained  during  these 
air  permeability  tests  are  shown  in  Table  3~B  of  Appendix  B. 

The  core  was  then  evacuated  for  2b  hours  and  at  the  end  of  this 
time  saturated  with  varsol.  Knowing  the  amount  of  varsol  injected,  the  porosity 
of  the  core  was  calculated.  Varsol  was  then  displaced  through  the  core  at  7 
different  rates  and  liquid  permeabilities  were  calculated.  The  data  for  these 
liquid  permeability  tests  are  shown  in  Table  of  Appendix  B.  Subsequently 

during  the  project ,  after  the  number  of  pressure  transducers  had  been  reduced 
to  3  from  4,  pure  varsol  was  again  displaced  through  the  core  at  5  additional 
rates.  At  the  time  of  these  latter  displacements  pure  propane  was  also  dis¬ 
placed  through  the  core  at  15  different  rates. 

Prior  to  the  initiation  of  an  actual  miscible  type  displacement  , 
the  clean  core  was  evacuated  for  2b  hours.  The  core  was  saturated  with  varsol 
and  a  check  on  the  core  porosity  was  obtained.  Approximately  5  pore  volumes 
of  crude  oil  were  then  injected  into  the  core  displacing  the  varsol.  If  some 
reservoir  fluid  other  than  crude  oil  were  to  be  used  then  the  initial  saturation 
with  varsol  was  not  carried  out.  When  the  core  was  saturated  with  reservoir 
fluid  all  lines  leading  to  the  transducers  and  the  Heise  gauge  were  purged  of 
gas  and  filled  with  liquid.  A  back  pressure  was  then  imposed  upon  the  system 
and  the  pressure  in  the  core,  transducers,  and  displacing  fluid  containers,  was 
gradually  built  up  to  the  desired  displacement  pressure.  This  back  pressure 
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was  generally  in  the  region  of  1600  p.s.i.g. 

All  flow  through  the  core  was  then  stopped  and  the  system  was 
allowed  to  come  to  equilibrium.  When  this  condition  was  attained  the  trans¬ 
ducers  were  zeroed  and  initial  pump  and  production  readings  were  taken.  The 
displacement  pomp  was  once  again  started  and  additional  reservoir  fluid  was 
injected  through  the  core.  This  was  continued  until  stable  pressure  differen¬ 
tials  from  the  Visicorder  were  obtained  along  the  entire  core  length.  When 
this  occurred  the  injection  line  was  transferred  from  the  reservoir  fluid 
container  to  the  displacing  fluid  container  with  out  stopping  the  displacement 
pump.  A  stop  watch  was  then  started  and  all  production  and  injection  data  were 
recorded  from  this  time  on.  Displacements  were  carried  on  as  long  as  were 
deemed  necessary  in  order  to  obtain  sufficient  data. 

At  the  end  of  the  displacement  the  injection  pump  was  stopped 
and  the  core  was  depressured  by  bleeding  off  the  back  pressure  regulator.  Any 
production  of  fluids  which  occurred  during  this  interval  was  noted  and  recorded. 
After  the  core  had  been  depressured  it  was  flushed  with  pentane  until  the  effluent 
appeared  to  be  pure  pentane.  During  this  flushing  operation  all  effluent  was 
collected  and  measured.  The  liquid  portion  of  the  effluent  was  then  allowed  to 
weather  for  a  few  days  and,  after  the  pentane  had  evaporated,  any  reservoir  oil 
which  had  been  produced  during  the  flushing  operation  was  measured  and  recorded. 

After  the  core  had  been  flushed  with  pentane,  pure  propane  was 
flowed  through  it  for  several  hours.  The  propane  was  then  followed  by  approxi¬ 
mately  150  cubic  feet  of  nitrogen.  The  core  was  then  dried  by  alternately 
flowing  air  through  it  and  then  imposing  a  vacuum  on  it  for  four  days.  An 
attempt  was  made  to  duplicate  the  entire  procedure,  from  initial  saturation  to 
cleaning,  during  each  displacement. 
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RESULTS  AND  DISCUSSION 


Theoretical  Results  and  Discussion  of  Critical  Calculations: 

A  comparison  of  the  measured  and  predicted  critical  temperatures 
and  pressures  of  twelve  different  hydrocarbon  mixtures  is  made  in  Table  1. 

The  methods  used  for  calculating  criticals  were  those  of  Organiek' sj  Kurata 
and  Katz  partially  modified  and  Eilerts.  The  calculated  criticals  of  all 
mixtures  with  the  exception  of ’,CH“P“843h  were  available  in  the  Literature 

The  standard  deviations  in  predicted  critical  temperatures  and 
pressures  were  calculated  and  appear  in  Table  2.  The  following  equation  was 
used  for  this  calculation: 


Standard  Deviation  = 


]  1  =  1 

where  N  =  number  of  mixtures. 


N 


T 


(observed  value  -  predicted  value)* 


N 


From  Table  2  it  would  appear  that  the  method  of  Organiek  is  superior  to  the 
others  in  predicting  critical  pressures.  For  predicting  critical  temperatures, 
Organiek' s  and  Eilerts'  correlations  compare  favorably  and  both  give  better 
results  than  the  Kurata-Katz  procedure.  These  conclusions,  of  course,  are  ma.de 
from  comparison  of  a  limited  cumber  of  samples. 

The  compositions  of  each  of  the  mixtures,  with  the  exception  of 
mixture  number  11,  are  tabulated  in  Appendix  A  along  with  sample  calculations 
of  each  method  of  predicting  criticals.  All  mixtures  contained  appreciable 
amounts  of  methane,  ranging  from  64$  in  mixture  C-l  to  82$  in  mixtures  B-l 
and  C“2.  The  amount  of  Cy+  components  in  each  mixture  was  relatively  small j 
from  2$  in  mixture  C-2  to  12$  in  several  of  the  other  mixtures.  Two  of  the 
mixtures,  CH-P~843  and  E,  contained  nitrogen  as  well  as  carbon  dioxide  while 
the  rest  contained  C0.p  only.  From  Table  1  it  can  be  seen  that  the  deviations 
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■were  large  for  all  methods  when  the  mixture  contained  nitrogen.  All  the 
mixtures  studied  could  be  classified  as  gas  condensates  because  of  the  small 
amount  of  liquid  present  at  equilibrium  conditions.  The  mole  per  cent  of 

the  mixture  in  the  vapor  phase  at  equilibrium  conditions  at  1^.7  p.s.i.a. 

o 

and  75  F  varied  from  85 $>  for  mixture  D  to  98 1°  fo r  mixture  C-2. 

Examining  each  method  of  predicting  critical s  individually,  a  few 
general  conclusions  can  be  drawn.  Organick's  correlation  appears  to  be  quite 
good  for  those  mixtures  with  high  methane  concentrations  and  little  or  no  acid 
gas  components.  However,  if  the  concentration  of  the  lighter  components 
becomes  too  large,  resulting  in  a  low  weight  average  equivalent  molecular  weight, 
Organick’s  correlation  becomes  inaccurate.  This  is  due  to  the  fact  that  the 
lines  on  the  correlation  grid  change  slope  rapidly  in  this  region  at  low  weight 
average  equivalent  molecular  weights.  The  accuracy  of  the  Organiek  correlation 
can  also  be  questioned  when  one  is  working  in  the  central  regions  of  the  grid 
where  all  lines  have  been  drawn  in  by  interpolation. 

The  calculations  involved  in  the  Organiek  prediction  are  relatively 
simple  providing  an  accurate  A.S.T.M.  distillation  analysis  of  the  Cl,  portion 
of  the  liquid  phase  is  available.  Without  this  distillation  analysis,  the 
accuracy  of  the  calculation  suffers. 

The  prediction  of  criticals  using  Eilert^  correlations  involve  the 
most  complex  calculations  of  any  of  the  three  methods  studied.  Not  only  is  an 
accurate  U.S.B.M.  distillation  of  the  Cy+  portion  of  the  mixture  required  but 
also  the  vapor-liquid  equilibrium  ratios  of  the  mixture  must  be  calculated  at  a 
standard  state.  This  of  course  entails  a  trial  and  error  procedure.  The  pre¬ 
diction  of  critical  temperatures  by  the  Eilerts  method  compares  quite  favorably 
with  results  obtained  using  Organick's  procedure.  This  is  because  both 


. 
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procedures  use  the  normal  absolute  boiling  point  of  each  component  as  a  weighting 
factor  when  calculating  critical  temperature. 

Prediction  of  critical  pressure  by  the  Eilerts  method  is  not  nearly 
as  accurate  as  prediction  of  critical  temperature  especially  if  use  is  made  of 
the  original  charts  published  in  Monograph  10  of  the  U.S.  Bureau  of  Mines.  While 
Eilerts  used  all  of  the  mixtures  presented  in  Appendix  A  to  develops  his  correla¬ 
tion,  he  balanced  positive  and  negative  deviations  when  drawing  the  final  charts. 
Consequently,  use  of  these  charts  can  sometimes  lead  to  a  considerable  error  in 

the  answer.  This  is  particularly  true  when  the  mole  fraction  of  the  mixture 

,  o 

in  the  vapor  phase  at  14.7  p.s.i.a.  and  75  F  approaches  1.0.  In  this  region  the 
slopes  of  the  lines  on  the  correlation  chart  approach  infinity  and  it  is  difficult 
to  obtain  a  satisfactory  answer.  This  problem  is  also  encountered  when  using 
the  critical  temperature  correlation  chart.  In  the  region  where  the  mole 
fraction  in  the  vapor  phase  is  less  than  0.5»  all  lines  have  been  extrapolated 
on  both  charts  and  the  accuracy  of  the  correlation  in  this  region  is  questionable. 

As  was  the  case  with  Organick's  method,  the  Eilerts  correlation 
gives  best  results  when  the  mole  fraction  of  the  mixture  in  the  vapor  phase  is 
below  0.90  and  when  the  mixture  contains  little  or  no  nitrogen.  The  secondary 
disadvantage  of  the  Eilerts  procedure  is  the  laborious  calculations  which  have 
to  he  carried  out  in  order  to  make  use  of  the  correlation.  It  is  doubtful  if 
the  accuracy  of  the  results  obtained  justify  the  work  involved. 

The  third  method  of  predicting  criticals  which  was  investigated 
wa.s  a  partially  modified  version  of  the  Kurata-Katz  procedure.  This,  as  fa.r  as 
can  be  determined,  was  the  method  used  by  Beriham,  Dowden  and  Kunzman  when 
they  developed  their  correlation  for  predicting  miscibility.  In  this  procedure 
use  is  made  of  the  M*  correction  developed  by  Davis,  Bertuzzi.,  Gore  and  Kurata 
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to  account  for  the  presence  of  nitrogen  in  the  mixture  •when  calculating 
critical  pressure.  However,  critical  temperatures  predicted  by  Benham, 

Dowden  and  Kunzman  were  obtained  using  the  original  Ku rata -Katz  procedure. 

The  predicted  criticals  shown  in  Table  1  under  the  heading  Ku rata -Katz  were 
therefore  calculated  in  the  same  manner. 

The  predicted  criticals  shown  in  Table  1  and  the  standard  deviations 
shown  in  Table  2  indicate  that  the  Kurata-Katz  method  was  the  least  accurate 
procedure  of  the  three  investigated.  Also  there  did  not  seem  to  be  any  trend 
to  the  Kurata-Katz  method  with  regard  to  mole  fraction  of  vapor  present  or 
mole  fraction  of  nitrogen  present.  One  drawback  to  predicting  critical 
pressures  with  this  method  is  the  fact  that  the  correlation  graph  plots  criti¬ 
cal  pressure  as  a  function  of  molal  average  critical  pressure,  molal  average 
critical  temperature,  true  calculated  critical  temperature  and  molal  average 
molecular  weight  of  the  mixture.  Thus  if  the  true  calculated  critical  tempera¬ 
ture  is  in  error  then  the  predicted  critical  pressure  obtained  from  the  correla¬ 
tion  chart  will  probably  also  be  in  error.  In  addition  to  this  the  molal  average 
molecular  weight  can  be  inaccurate  unless  sufficient  information  about  the  C^+ 
portion  of  the  mixture  is  available. 

The  advantage  of  the  Kurata-Katz  procedure  is  its  relative  simpli¬ 
city  and  the  comparatively  easy  calculations  which  are  involved.  When  a  dis¬ 
tillation  analysis  of  the  portion  of  the  mixture  is  unavailable  then  this 
method  or  its  modified  form  is  definitely  the  best  procedure  to  use  when  calcu¬ 
lating  critical  conditions.  The  limitations  which  are  involved  when  using 
this  method  must  be  kept  in  mind  however. 

The  accuracy  of  all  the  procedures  available  for  calculating 
critical  conditions  is  necessarily  dependent  upon  the  accuracy  of  the  mixture 
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analysis.  It  is  therefore  imperative  that  as  detailed  an  analysis  as  possible 
be  obtained,  if  calculations  of  the  critical  temperature  and  pressure  of  a 
naturally  occurring  hydrocarbon  mixture  are  to  be  attempted. 


-  33  - 


Discussion  of  Experimental  Results; 

Initial  displacements  were  carried  out  on  the  Berea  sandstone 
core  to  determine  both  its  air  and  liquid  permeabilities  as  well  as  the 
porosity  of  the  core.  These  initial  displacements  also  served  to  familiarize 
the  investigators  with  the  experimental  equipment. 

Air  was  flowed  through  the  core  at  10  different  rates  and  the 
permeability  of  the  core  to  air  was  calculated  at  each  of  these  rates.  The 
results  of  these  calculations  are  shown  in  Table  3-B.  From  these  observations 
the  average  permeability  of  the  core  to  air  was  found  to  be  293.9  millidarcies. 
An  attempt  was  made  to  plot  air  permeability  against  the  inverse  of  the 
mean  pressure  in  order  to  illustrate  the  Klinkenberg  effect.  However ,  the 
wedely  scattered  data  which  were  obtained  for  the  low  flow  rates  used  made 
this  impractical. 

The  liquid  permeability  of  the  core  was  obtained  by  flowing 
varsol  through  it  at  various  rates  and  absolute  pressures.  The  data  resulting 
from  these  displacements  are  shown  in  Table  4-B.  In  all,  37  displacements 
were  carried  out  at  different  reservoir  conditions.  The  resulting  average 
liquid  permeability  of  the  core  was  calculated  to  be  266.4  millidarcies.  While 
there  was  some  variance  in.  the  experimental  data  obtained,  all  of  the  runs 
were  used  when  the  average  permeability  was  calculated.  The  most  pronounced 
deviations  were  observed  to  occur  at  low  absolute  reservoir  pressures  during 
low  flow  rates.  As  the  reservoir  pressure  was  Increased  during  these  low 
rates,  the  permeabilities  obtained  seemed  to  become  more  stable.  The 
tendency  for  the  permeability  to  vary  with  reservoir  pressure  was  not  as 
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marked  at  higher  displacement  rates.  Prom  these  results  it  would  appear 
that  liquid  permeabilities  should  be  obtained  during  relatively  high 
displacement  rates  (200  cc  per  hour  or  greater)  for  a  porous  media  of  the 
type  -under  investigation. 

Additional  displacements  of  pure  fluids  through  the  core  were 
carried  out  at  various  times  during  the  investigation.  To  preserve  continuity 
all  of  these  "pure  fluid  displacements"  will  be  discussed  at  this  time. 

Tables  5-B  to  8~B  inclusive  show  the  pressure  data  which  were  obtained  during 
displacements  of  varsol,  propane  and  pentane  through  the  core.  The  data  in 
Table  5“B  were  obtained  while  four  differential  pressure  transducers  were 
mounted  on  the  core.  The  data  in  the  other  three  tables  were  taken  using 
only  three  transducers.  The  data  from  each  of  the  Tables  were  plotted  and 
appears  as  Figures  1-B  to  4~B  inclusive. 

All  results  indicated  that  the  pressure  differentials  across 
the  entire  core  became  increasingly  non-uniform  as  the  displacing  rate 
increased.  That  is  to  say  adjacent  sections  of  the  core  were  not  subjected 
to  the  same  pressure  drops ,  especially  at  higher  flow  rates.  This  of  course 
would  be  natural  if  the  permeability  of  the  core  was  not  uniform  throughout 
its  entire  length.  When  varsol  and  propane  were  used  as  the  Injection  fluids, 
section  2,  or  the  middle  portion  of  the  core,  consistantly  was  subjected 
to  the  least  pressure  differential.  This  would  indicate  that  section  2  was 
slightly  more  permeable  than  the  end  sections  of  the  core. 

The  data  shown  for  the  displacement  of  pure  pentane  through  the 
core  were  obtained  at  a  later  date  than  the  data  for  either  the  varsol  or  pro¬ 
pane  displacements.  The  pentane  displacements  were  carried  out  approximately 
midway  through  the  Investigation,  prior  to  the  start  of  Displacement  No.  F. 
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After  completion  of  Displacement  No.  E  and  before  the  pentane  displacements 
took  place,  the  core  ms  reversed  in  the  flanges  so  that  the  original  outlet 
end  became  the  inlet  face.  The  reasons  for  this  reversal  will  be  discussed 
later.  The  pressure  differentials  obtained  while  flowing  pentane  in  this 
reversed  direction  are  shown  on  Figure  4-B.  On  this  plot  the  pressure 
differentials  across  sections  1  and  2  are  quite  similar  while  the  pressure 
drop  across  section  3  is  greater  than  for  either  of  the  first  two  sections. 

This  is  increasingly  apparent  as  the  flow  rates  increase.  This  phenomena 
could  be  due  to  damage  to  the  original  inlet  end  of  the  core  during 
Displacements  No.  D  and  No.  E. 

Five  miscible,  liquid-liquid  displacements  were  run  on  the  core 
during  the  investigations.  These  displacements  were  numbered  D,  E,  Fs  G, 
and  H.  The  production  data  for  each  of  them  is  contained  in  Tables  9~B  to 
13-B  inclusive.  Pressure  profile  data  was  obtained  on  four  of  the 
displacements,  D,  E,  F  and  H  and  is  recorded  in  Tables  l4-B  to  1‘7-B 
inclusive.  For  each  run,  with  the  exception  of  number  G,  plots  were  drawn 
of:  l)  Cumulative  Oil  Production  ($>  OIP)  versus 
Cumulative  Fluid  Injected  (pore  volumes), 

2)  Gas-Oil  Ratio  (cc/cc)  and  Oil  Production 

Rate  (ce/min)  versus  Cumulative  Fluid  Injected  (pore  volumes), 
and  3)  The  Pressure  Drop  per  Inch  at  any  time  versus  the  Original  Pressure 
Drop  per  Inch  (p. s. i.d./in. ). 

A  plot  of  this  last  relationship  is  not  included  for  Displacement  No.  G  due  to 
the  low  pressure  differentials  which  were  present  during  this  ran.  The  results 
of  these  miscible  displacements  will  be  discussed  in  the  chronological  order 


in  which  they  were  run. 
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The  initial  miscible  flood  performed  ms  Displacement  No.  D 

where  crude  oil  was  displaced  by  liquid  propane  at  a  rate  of  60  ec  per 

hour.  The  ratio  of  the  viscosity  of  the  crude  oil  to  the  viscosity  of  the 

liquid  propane  was  70.5.  During  this  flood,  a  recovery  to  initial  solvent 

breakthrough  of  approximately  30  per  cent  of  the  original  oil  in  place  was 

obtained.  The  recovery  of  oil  at  secondary  solvent  breakthrough  was  55 

per  cent  of  OIP  and  the  final  recovery  was  75  per  cent  of  OIP  after  injection 

of  approximately  1.7  pore  volumes  of  liquid  propane.  These  recoveries;,  while 

seemingly  low  for  a  miscible  type  drive,  are  in  general  agreement  with  results 

(k  24) 

obtained  by  other  investigators*  *  .  The  poor  recoveries  appear  to  be 

caused  mainly  by  the  adverse  viscosity  ratio  encountered  during  the  displacement 
of  a  relatively  viscous  crude  oil  by  liquid  propane. 

A  better  idea  of  the  displacement  process  can  be  obtained  from 
the  plot  of  GOR  and  oil  production  rate  versus  pore  volumes  of  fluid  injected. 

Up  to  injection  of  approximately  0.28  pore  volumes,  the  effluent  produced 
to  the  separator  was  a  single  liquid  phase.  At  this  point  initial  solvent 
breakthrough  occurred  and  recovery  of  two  separate  phases  began.  After 
approximately  0.68  pore  volumes  of  propane  had  been  injected,  a  suggested 
second  solvent  breakthrough  occurred  as  indicated  by  a  rapid  rise  in 
GOR  and  a  corresponding  decline  in  production  rate.  The  first  premature 
arrival  of  propane  may  he  attributed  to  the  effects  of  viscous  fingering. 

From  this  first  solvent  arrival  until  the  secondary  breakthrough,  it  is 
suggested  that  the  propane  and  oil  are  intimately  mixed  and  completely  miscible 
such  that  the  fluid  flowing  is  an  oil  enriched  with  injected  propane.  There 
is  some  evidence  from  the  results  that  the  mechanism  after  the  secondary 
breakthrough  is  similar  to  the  subordinate  stripping  action  of  a  non-equilibrium 
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two-phase  gas  drive  displacement  thus  accounting  for  the  high  residual  oil 
saturation  at  the  termination  of  the  test. 

The  pressure  profile  plot  for  this  displacement  also  indicates 
the  presence  of  viscous  fingers.  The  profile  along  the  core  is  initially 
stabilized  and  is  a  straight  line  at  unit  value.  With  the  injection  of 
displacing  fluid,  a  zone  of  very  low  pressure  drop  advanced  into  the  core. 
This  zone  of  low  pressure  drop  reached  the  downstream  6  inches  of  the  core 
after  only  0.139  pore  volumes  of  fluid  had  been  injected.  The  pressure 
drop,  of  course,  was  primarilly  due  to  the  highly  viscous  oil  being  displaced 
by  a  fluid  of  much  lower  viscosity. 

The  second  miscible  flood  undertaken  was  Displacement  No.  E 
where  commercial  varsol  was  displaced  by  liquid  propane  at  a  rate  of  200 
cc  per  hour.  The  ratio  of  the  viscosity  of  the  varsol  to  the  viscosity 
of  the  propane  in  this  case  was  8.6.  The  recovery  of  original  varsol  was 
45  per  cent  of  "initial  oil  in  place"  when  solvent  breakthrough  occurred 
during  this  flood.  At  the  end  of  the  test  approximately  1.8  pore  volumes 
of  propane  had  been  injected  into  the  core  and  95  per  cent  of  the  "initial 
oil  in  place"  had  been  recovered.  (During  subsequent  discussions  the  term 
"initial  oil  in  place"  shall,  refer  to  whatever  fluid  the  core  was  saturated 
with  before  a  displacement  was  carried  out).  During  blowdown  and  cleanup 
operations,  after  completion  of  the  run,  all  remaining  varsol  in  the  core  was 
recovered. 

Neither  the  recovery  data  nor  the  GOR  and  oil  production  rate 
data  indicated  the  occurance  of  a  secondary  solvent  breakthrough  during 
this  displacement.  Initial  solvent  breakthrough  occurred  after  0.4  pore 
volumes  of  propane  had  been  injected  into  the  core.  From  this  point  until 
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the  end  of  the  test  the  GOR  steadily  increased  -while  the  oil  production  rate 
decreased  gradually.  These  results  indicate  that  viscous  fingering  was  not 
such  a  prominant  factor  during  this  displacement  as  it  was  during  Displacement 
No.  D.  This  is  also  borne  out  by  the  fact  that  the  pressure  profile  data 
shows  that  the  outlet  12  inches  of  the  core  did  not  exhibit  a  drop  is  pressure 
differential  until  0.26l  pore  volumes  of  propane  had  been  injected. 

After  completion  of  Displacement  Ho.  E,  an  attempt  was  made  to 
rerun  the  previous  test  involving  crude  oil  and  propane .  Difficulties  were 
encountered  however  when  the  core  was  being  saturated  with  crude  oil.  During 
this  saturation  it  was  observed  that  practically  all  of  the  pressure  drop 
across  the  core,  when  crude  oil  was  flowing  at  60  cc  per  hour,  was  confined 
to  the  inlet  12  inches  of  the  core.  The  core  was  baekf lushed  with  pentane 
and  cleaned  and  dried  before  another  attempt  was  made  to  saturate  it  with 
crude  oil.  Again  the  results  were  unfavorable  with  the  definite  indication 
of  plugging  taking  place  at  the  inlet  end.  The  core  was  again  baekf lushed  and 
cleaned  and  was  then  dismantled  from  its  flanges.  Upon  examining  the  core  it 
was  observed  that  the  rubber  end  gaskets  had  swelled  until  they  were  partially 
covering  the  face  of  the  core.  It  was  also  noted  that  the  inlet  end  of  the  core 
was  very  heavily  oil-stained  particularly  on  its  bottom  side.  This 
indicated  the  possibility  that  gravity  segregation  took  place  during  the  flow 
of  the  heavy  crude  oil.  It  was  decided  to  thoroughly  clean  the  core  with 
pentane,  replace  the  original  rubber  gaskets  with  new  gaskets  and  remount 
the  core  in  the  flanges  in  a  reversed  position.  This  was  done  with  the 
result  that  the  original  partially-damaged,  inlet  end  now  became  the  outlet 
end.  From  this  time  on  it  was  noticed  that  the  pressure  differential  across 
section  3  i.  e.  the  outlet  section,  was  slightly  higher  than  it  was  across  either 


of  the  other  two  sections. 
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The  next  miscible  flood  undertaken  was  Displacement  No.  F 
which,  along  with  Displacement  No.  G,  represent  the  portion  of  the  investi¬ 
gation  carried  out  at  low  viscosity  ratios.  Displacements  No.  F  and  No.  G 
both  entailed  the  displacement  of  commercial  pentane  by  liquid  propane  at 
rates  of  400  and  200  cc  per  hour  respectively.  The  ratio  of  the  viscosity  of 
pentane  to  the  viscosity  of  propane  in  these  cases  was  2.13.  The  recoveries 
of  original  pentane  obtained  at  solvent  breakthrough  were  JO.O  per  cent  of 
OIP  at  the  faster  rate  and  65.0  per  cent  of  OIP  at  the  slower  rate.  Final 
liquid  pentane  recoveries  of  79  and  jQ  per  cent  of  OIP  were  obtained  after 
only  0.9  pore  volumes  of  propane  had  been  injected.  From  this  point  until 
completion  of  the  test,  no  additional  liquid  pentane  was  recovered  during 
either  displacement.  Material  balance  calculations  were  carried  out  on 
both  displacements  and  they  indicated  that  approximately  20$  of  the  original 
pentane  in  the  core  was  vaporized  in  the  separator  and  was  measured  in  the 
effluent  as  a  gas  along  with  the  propane.  On  the  basis  of  these  material 
balance  calculations  and  the  fact  that  during  the  cleanup  of  the  core  no 
additional  liquid  pentane  was  recovered,  it  can  he  assumed  that  100  per  cent 
of  the  original  pentane  in  place  was  recovered  during  both  displacements. 

Up  to  the  time  of  initial  solvent  breakthrough  it  is  thought  that  very  little 
pentane  was  flashed  into  the  vapor  phase.  However,  after  arrival  of  the 
miscible  liquid  to  the  separator,  it  is  believed  that  portions  of  the 
liquid  pentane  started  to  vaporize  as  the  propane  flashed  from  liquid  to  gas. 

Again,  as  in  the  case  of  propane  displacing  varsol,  neither  the 
recovery  data  nor  the  GOR  and  oil  production  data  indicated  the  occur ranee 
of  a  secondary  solvent  breakthrough.  Initial  solvent  breakthrough  occurred  at 
the  high  rate  after  0.68  pore  volumes  of  propane  had  been  injected.  At  the 
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low  rate,  initial  solvent  breakthrough  was  observed  after  0.64  pore  volumes 
of  propane  had  been  injected.  It  is  possible  that  at  the  lower  rate, 
diffusional  effects  at  the  flood  front  between  the  propane  and  pentane  could 
account  for  an  earlier  breakthrough  and  resultant  smaller  recovery  of 
initial  oil  in  place.  There  was  no  evidence  of  viscous  fingering  taking  place 
during  either  Displacements  Wo.  F  or  Wo.  G.  This,  of  course,  is  to  be 
expected  when  the  viscosity  ratio  of  the  displaced  and  displacing  fluids  is 
relatively  low.  The  pressure  profile  data  for  Displacement  Wo.  F  confirmed 
this  as  the  outlet  12  inches  of  the  core  did  not  show  a  drop  in  pressure 
differential  until  0.556  pore  volumes  of  propane  had  been  injected  into  the 
core. 

The  final  miscible  flood  performed  was  Displacement  Wc.  H 
where  a  heating  oil,  Esso  P-40,  was  displaced  by  liquid  propane  at  a  rate 
of  200  cc  per  hour.  The  ratio  of  the  viscosity  of  the  heating  oil  to  the 
viscosity  of  the  liquid  propane  was  14.2.  During  this  displacement,  a 
recovery  to  initial  solvent  breakthrough  of  appro xlmately  35  per  cent  of 
the  original  oil  in  place  was  obtained.  As  was  the  case  with  Displacement 
Wo.  D  where  crude  oil  was  displaced  by  propane,  it  is  thought  that  a  secondary 
solvent  breakthrough  was  detected.  The  recovery  of  oil  at  this  secondary 
breakthrough  was  approximately  57  per  cent  of  0IP.  The  final  recovery  obtained 
during  this  flood  was  76  per  cent  of  0IP  after  injection  of  1.5  pore  volumes 
of  liquid  propane. 

The  similarity  between  this  heating  oil  displacement  and  the 
crude  oil  displacement  Is  very  evident  even  though  there  is  a  five  fold 
difference  in  viscosity  ratios  between  the  two  floods.  During  Displacement 
Wo.  H,  the  heating  oil  test,  the  effluent  produced  to  the  separator  was 
a  single  liquid  phase  up  to  injection  of  0.33  pore  volumes  of  propane.  From 
this  point  on  recovery  of  two  separate  phases  was  observed.  After  injection 


of  0.71  pore  volumes  of  propane,  a  rapid  rise  in  GOR  occurred  accompanied 
by  a  corresponding  decrease  in  oil  production  rate.  This  suggested  the 
occurrance  of  a  secondary  solvent  breakthrough  similar  to  the  one  experienced 
during  Displacement  No.  D.  Again  it  is  suggested  that  the  flood  exhibited 
miscible -type  characteristics  up  to  the  time  of  secondary  solvent  breakthrough, 
but  after  this  time  the  displacement  resembled  the  subordinate  stripping 
action  of  a  non-equilibrium,  two -phase,  gas  drive.  The  pressure  profile  data 
for  this  displacement  also  suggested  the  presence  of  viscous  fingers.  A 
decrease  in  pressure  differential  across  the  outlet  12  inches  of  the  core  was 
observed  after  injection  of  only  0.216  pore  volumes  of  liquid  propane. 

A  summary  of  all  recovery  data  obtained  during  the  investigation 
appears  in  Table  3. 


TABLE  3 

RECOVERIES  OF  RESERVOIR  LIQUID 


Rec.  of  Res. 


Displacement 

Displacement 

Viscosity  Ratio 

Rec.  of  Res. 

Liquid  After 

Number 

Rate 

^oil/  M- solvent 

Liquid  at  B.T. 

1.0  FV  Dis¬ 
placed 

(cc/hr)  ($  OIP) 


No.  D 

60 

70.5 

30 

65 

No.  E 

200 

8.6 

45 

83 

No.  F 

400 

2.13 

TO 

79.5  + 

No.  G 

200 

2.13 

65 

78  + 

No.  H 

200 

14.2 

35 

67 

The  Cumulative  Oil  Recovery  to  Initial  Solvent  Breakthrough  was  plotted  versus 
the  Viscosity  Ratios  of  the  oils  and  solvent  and  appears  as  Figure  6.  Super- 

(24) 

imposed  upon  this  plot  are  the  results  obtained  by  Offeringa  and  va,n  der  Peel' 
during  their  miscible  and  immiscible  tests  on  scaled  models.  As  can  be  seen 
there  is  relatively  good  agreement  between  their  results  and  the  data  obtained 
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during  this  investigation.  The  greatest  degree  of  divergency  occurs  in 
the  region  of  high  viscosity  ratios. 

It  is  of  interest  to  note  that  in  the  region  of  viscosity 
ratios  greater  than  10.0,  the  change  in  recovery  of  oil  at  solvent  break 
through  is  relatively  small  for  any  marked  change  in  viscosity  ratio.  In. 
this  region  of  adverse  viscosity  ratios,  it  is  evident  that  the  early 
solvent  breakthroughs  are  a  result  of  the  formation  of  viscous  fingers 
during  the  displacement.  In  the  region  of  viscosity  ratios  between  1.0 
and  10.0  however,  recovery  of  oil  to  solvent  breakthrough  changes  appreciably 
with  a  change  in  viscosity  ratio.  This,  of  course,  indicates  a  constantly 
decreasing  effect  on  recovery  due  to  viscous  fingering  as  the  viscosity  ratio 
decreases  to  unity.  Therefore  when  the  ratio  of  the  viscosity  of  the 
reservoir  oil  to  the  viscosity  of  the  displacing  fluid  is  of  a  magnitude 
less  than  10.0  it  is  believed  that  the  displacement  mechanism  which  takes 
place  is  primarily  of  a  " plunger”  type,  where  a  graded  transition  zone  is 
formed  between  the  oil  and  the  displacing  fluid.  However  at  high  viscosity 
ratios  the  displacing  fluid  travels  through  the  reservoir  quite  rapidly 
developing  a  wide  transition  zone  with  many  fingers.  Thus  at  high  viscosity- 
ratios,  any  attempt  to  correlate  some  "length  of  mixing  zone"  is  meaningless. 

These  conclusions  are  somewhat  verified  by  the  data  presented 
in  Table  4.  Here  the  "rate  of  advance  of  the  leading  edge  of  the  misting  zone" 
is  tabulated  for  the  different  displacements  performed.  These  data  were 
obtained  from  Visieorder  charts.  Comparing  the  tests  which  were  carried  out 
at  the  same  displacement  rate,  i.e.  Displacements  No.  E,  G  and  H,  a  variance 
in  the  rate  of  advancement  of  the  mixing  zone  can  be  seen,  especially  across 
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the  first  two  sections  of  the  core.  During  Displacement  No.  G  when  the 
viscosity  ratio  was  2.13,  the  mixing  zone  stabilized  quickly  and  travelled 
through  sections  one  and  two  at  exactly  the  same  rate.  However,  as  the 
viscosity  ratio  increased,  the  mixing  zone  spread  initially  quite  quickly 
through  section  one  and  then  progressively  slowed  down  as  it  advanced 
toward  the  outlet  face.  This  initial  rapid  advancement  became  more  pro¬ 
nounced  as  the  viscosity  ratio  increased  from  8.6  to  14.2.  These  limited 

( 20) 

observations  are  in  agreement  with  views  expressed  by  other  authors' 
when  they  investigated  this  phenomena  of  mixing  zone  growth.  In  order  to 
rigorously  examine  the  problem  of  mixing  zone  growth  it  is  suggested  that  a 
longer  porous  system  than  the  one  used  during  this  investigation  is  needed. 

A  visual  examination  of  the  interaction  between  liquid  propane 

and  Esso  P-40  heating  oil  was  made  in  a  laboratory  high  pressure  cell. 

o 

These  observations  were  carried  out  at  75  F  and  pressures  ranging  from  .100 
p.s.i.g.  to  2000  p.s.i.g.  Under  static  conditions,  an  interface  between  the 
propane  and  the  heating  oil  was  observed  to  be  present  throughout  the  entire 
range  of  pressures  investigated.  The  transition  zone  between  the  propane 
and  the  heating  oil  was  very  narrow  with  little  apparent  mixing  taking  place 
between  the  two  liquids.  However,  upon  agitation,  the  two  liquids  became  a 
single  fluid  even  at  the  lowest  pressure  observed.  No  difficulty  should  be 
encountered  therefore,  in  getting  two  fluids  of  this  type  to  mix  raiscibly 
while  travelling  through  a  porous  media. 

The  data  which  has  been  obtained  during  this  investigation  has 
resulted  in  the  formulation  of  some  general  conclusions  regarding  the  mechan¬ 
isms  which  occur  during  a  miscible  type  displacement.  It  is  felt  that  the 
efficiency  of  any  laboratory  miscible  displacement  (compared  to  that  predicted 
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theoretically)  will  be  primarily  influenced  by  the  sweep  efficiency  within 
the  core.,  This*  in  turn,  is  controlled  by  such  factors  as  unfavorable 
mobility  or  viscosity  ratio  effects,  reservoir  inhomogeneities  and  gravitational 
segregation  effects.  Thus  attempts  to  verify  the  predicted  miscibility 
conditions  of  a  system  by  laboratory  tests  will  necessitate  the  suppress- 
ment  or  elimination  of  all  causes  of  poor  sweep  efficiency.  In  this  regard, 
gravity  effects  may  be  eliminated  by  using  vertical  systems  while  the 
development  of  fingering  can  perhaps  be  controlled  by  carrying  out  displace¬ 
ments  at  or  below  some  critical  velocity. 
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Discussion  of  Experimental  Apparatus; 

A  brief  discussion  •will  be  presented  at  this  time  outlining  the 
advantages  and  disadvantages  of  the  experimental  apparatus  which  was  designed 
for  this  investigation.  The  most  delicate  components  of  the  apparatus  were 
the  differential  pressure  transducers.  Originally  four  in  number s  they  were 
mounted  in  series  on  a  rack  below  the  core.  Stainless  steel  tubings  one 
eighth  inch  in  diameter,  was  used  to  connect  the  pressure  taps  on  the  core 
to  the  high  and  low  sides  of  each  transducer.  Each  pressure  transmitting 
line  was  valved  once  between  the  core  and  the  transducer.  After  considerable 
difficulty  had  been  encountered,  it  was  realized  that  a  line  connected  the 
ends  of  each  transducer  was  necessary.  The  purpose  of  this  equalizing 
connection  was  to  maintain  the  same  pressure  on  each  side  of  a  transducer 
during  the  time  the  entire  system  was  being  pressured  up  or  depressured. 

This  type  of  arrangement  prevented  any  sudden,  unsafe  pressure  drops  being 
imposed  across  a  transducer  particularly  during  the  period  when  the  core  was 
being  depressured.  When  the  number  of  transducers  in  use  was  reduced  to  three, 
the  two  end  transducers  were  remounted  in  a  manner  just  described  with  satis¬ 
factory  results. 

The  extreme  care  which  is  required  when  calibrating  and  using 
this  type  of  pressure  sensing  device  cannot  be  overstressed.  If  additional 
transducers  were  added  to  the  system,  the  entire  bank  of  instruments  will 
probably  have  to  be  mounted  in  a  location  remote  from  the  core  itself.  This 
should  present  no  difficulty  or  resulting  loss  in  accuracy.  The  full  scale 
deflection  of  the  diaphram  in  the  transducer  requires  the  transmission  of 
only  about  0.1  cc  of  liquid  with  the  result  that  pressure  changes  in  the 
core  are  detected  almost  instantaneously  providing  the  pressure  Hines  axe 
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filled  -with  a  non-compressible  liquid. 

The  operation  of  the  Visicorder  presented  little  difficulty 
once  the  investigators  became  familiar  with  the  instrument.  The  galvanometers 
were  zeroed  prior  to  the  start  of  any  test  and  then  rechecked  after  the  core 
and  transducers  were  pressured  up  to  the  desired  reservoir  conditions.  Care 
was  taken  during  this  operation  that  the  voltage  supply  to  the  transducers 
was  constant  and  stabilized.  Occasionally  the  traces,  corresponding  to  the 
two  upstream  galvanometers  were  noticed  to  vibrate  rapidly  during  the  preparation 
for  a  displacement.  It  was  subsequently  discovered  that  this  vibration  was 
caused  by  the  transmission  of  harmonics  which  were  present  in  the  relatively 
long  feed  line  from  the  injection  pump  to  the  core.  These  vibrations  could 
be  damped  out  immediately  by  attaching  the  inlet  feed  line  securely  'bo  the 
frame  of  the  pump.  Adjustment  must  be  made  to  the  Visicorder  every  time  a 
different  paper  speed  is  used  to  ensure  that  the  trace  intensification  of  the 
instrument  remains  uniform.  This  relatively  easy  adjustment  will  result  in 
clearer  and  sharper  traces  being  obtained.  All  records  from  the  Visicorder 
may  be  permanized  to  prevent  fading  by  a  relatively  simple  fixing  process. 
Although  this  was  not  done  during  the  course  of  this  investigation  it  is 
recommended  as  standard  procedure  for  any  subsequent  work. 

With  regard  to  the  remainder  of  the  test  apparatus,  three  main 
areas  of  difficulty  were  found.  Firstly,  the  accuracy  and  efficiency  of 
the  gas  measuring  and  venting  system  left  room  for  Improvement.  It  was 
discovered  early  in  the  investigation  that  the  burette-type  measuring  devices 
were  too  small  to  handle  the  large  quantity  of  propane  which  was  produced 
during  a  test.  Consequently  a  wet  test  meter  was  used  to  measure  all  gas  flows. 
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The  chief  disadvantage  of  this  procedure  was  the  inherent  inaccuracy  of  a 
wet  test  meter  during  low  flow  rates  and  its  inability  to  instantaneously 
detect  an  initial  gas  breakthrough.  Gas  breakthrough  therefore,  had  to  be 
observed  visually  in  the  separator  column  and  was  recorded  at  that  time  when 
the  first  bubble  of  gas  appeared  at  the  outlet  flowline  into  the  separator. 

The  system  provided  for  venting  the  produced  gas  was  also  inadequate,  mainly 
due  to  its  remote  location  from  the  test  apparatus.  Because  no  suction  can  be 
applied  to  the  discharge  line  from  the  wet  test  meter,  a  large  diameter  vent 
line  is  required  in  order  that  proper  dissipation  of  the  flammable  gas  can  be 
achieved. 

A  second  source  of  difficulty  during  the  tests  was  the  Grove 
Back -pressure  Regulator.  Designed  to  hold  a  constant  down  stream  pressure  on 
the  system,  it  was  found  that  after  breakthrough  of  propane  the  regulator  behaved 
very  erratically.  This  was  partly  due  of  course,  to  the  extreme  cooling  effect 
which  took  place  across  the  regulator  when  the  propane  flashed  into  a  vapor. 

When  freezing  did  occur  in  the  regulator,  it  was  easily  recognized  hy  a 
decrease  to  zero  in  production  and  a  corresponding  increase  in  pressure  in  the 
core  if  the  injection  pump  was  still  running. 

During  the  running  of  a  test,  both  the  upstream  and  downstream 
pressure  on  the  core  had  to  be  obtained  from  a  single  Heise  gauge.  This 
operation  entailed  the  constant  opening  and  closing  of  valves  and  probably 
resulted  in  some  inaccuracies  in  the  pressure  readings.  A  response  on  the 
Visicorder  chart  would  also  occur  whenever  a  valve  close  to  the  downstream 
end  of  the  core  was  adjusted.  For  these  reasons,  the  employment  of  two 
separate  Heise  gauges  on  the  system  could  probably  be  justified.  As  was 
pointed  out  earlier  in  the  discussion,  plugging  of  the  upstream  end  of  the 


. 


'  v'jfc  •  jart  e .lit  n.  »  i  l  ?  ixiW 

> 


-  49  - 


core  was  detected  after  the  repeated  use  of  crude  oil  as  a  reservoir  fluid. 
This  condition  could  he  avoided  by  mounting  a  filter  apparatus  between  the 
upstream  end  of  the  core  and  the  injection  pump. 
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SUMMARY  AM)  CONCLUSIONS 

It  can  be  concluded  from  the  experimental  results  that  the  most 
important  factor  influencing  the  results  of  a  laboratory  miscible  displacement 
test  is  the  sweep  efficiency  within  the  core.  The  sweep  efficiency,  in  turn,  is 
primarily  affected  by  adverse  mobility  or  viscosity  ratios,  reservoir  inhomo¬ 
geneities  and  gravity  segregation.  Of  these  three,  it  would  appear  that  the 
viscosity  ratio  was  the  dominant  factor  during  the  miscible  displacement  of 
reservoir  fluids  from  the  relatively  homogeneous  sandstone  used  during  this 
investigation.  Above  viscosity  ratios  of  10.0  the  presence  of  viscous 
fingers  was  indicated  both  by  early  solvent  breakthrough  and  by  rapidly 
decreasing  pressure  differentials.  At  viscosity  ratios  less  that  10.0  the 
limited  results  obtained  indicated  that  the  mechanism  of  frontal  advance  more 
closely  resembled  the  piston-like  displacement  visualized  by  many  investigators. 
However,  even  at  high  viscosity  ratios  where  the  presence  of  viscous  fingers 
was  suspected,  the  recovery  to  breakthrough  is  similar  in  magnitude  to  that 
obtained  by  water  floods  for  the  same  viscosity  ratio.  The  subordinate 
production  of  a  miscible  flood  is  much  greater  however  due  to  the  more 
efficient  stripping  action  which  takes  place.  It  is  possible  that  at  low 
viscosity  ratios  and  low  displacements  rates,  the  effect  of  diffusional 
mixing  at  the  flood  front  will  influence  the  recovery  of  oil  to  solvent 
breakthrough . 

The  theoretical  investigation  reveals  that  while  the  method  of 
predicting  criticals  used  by  Benham  et.  al.  during  development  of  their  correlation 
is  not  the  most  accurate  available,  it  Is  probably  reliable  enough  for  use  in 
an  approximate  correlation  such  as  theirs.  At  any  rate,  it  is  believed  that 
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the  controlling  factor  in  predicting  the  critical  conditions  for  a  certain 
fluid  system  in  a  porous  media  is  not  the  method  used  to  calculate  criticals 
but  is,  in  fact,  the  method  used  to  determine  the  composition  of  fluid  upon 
■which  the  calculations  are  made.  In  this  regard,  it  is  suggested  that  the 
equipment  designed  and  used  during  this  investigation  could  be  used  to  deter¬ 
mine  these  compositions,  particularly  during  miscible  liquid-liquid  dis¬ 
placements. 

In  the  case  of  liquid-liquid  displacements,  the  pressure 
differential  present  in  any  section  of  the  core  under  steady  state,  laminar 
flow  conditions,  is  directly  related  to  the  viscosity  of  the  fluid  which 
is  flowing  at  that  particular  time  and  in  that  particular  section.  Therefore, 
knowing  the  pressure  differential  across  any  portion  of  the  core,  it  would  be 
possible  to  calculate  the  viscosity  of  the  fluid  which  is  flowing.  Then  by 
knowing  the  viscosities  of  various  mixtures  of  the  displaced  and  displacing 
fluid,  the  composition  of  the  mixture  in  the  core  could  be  obtained.  This 
composition  profile  could  be  correlated  directly  to  the  pressure  differential 
which  would  occur  in  the  core  during  the  flow  of  one  of  these  mixtures. 

Perhaps  the  simplest  method  of  doing  this  would  be  to  flow  first,  pure  fluid 
through  the  core,  observing  the  pressure  drops,  and  then  flow  various  known 
mixtures  through  the  core  and  compare  pressure  differentials.  The  viscosities 
of  these  known  mixtures  could  be  calculated  from  the  pressure  differentials 
and  checked  in  a  laboratory  viscometer.  This  procedure  would  also  provide  a 
comparison  between  what  occurs  in  an  idealized  laboratory  cell  versus  the 
actual  behavior  in  a  porous  media.  The  pressure  transducers  and  Visicorder 
equipment  would  be  useful  in  many  other  phases  of  reservoir  studies.  Such 
things  as  water  floods,  alcohol  displacements,  gas  displacements  or  any  two- 
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phase  operation  could  employ  the  use  of  equipment  such  as  this. 

While  the  original  purpose  of  this  investigation  -was  to  check 
the  accuracy  of  theoretical  methods  for  predicting  critical  conditions ,  it 
was  soon  discovered  that  the  confimation  of  predicted  miscibility  pressures  and 
temperature  by  laboratory  testing  procedures  necessitated  the  elimination  or 
suppressment  of  all  cause  of  poor  pore  displacement  efficiency.  The  major 
portion  of  the  work  involved  in  this  investigation  was  therefore  directed 
toward  this  endeavour. 

Restating  the  Conclusions; 

1)  The  method  of  predicting  eritie&ls  conditions  developed  hy 
Benham,  Dowden  and  Kunzman  should  he  used  as  a  first  approximation  only. 

2)  The  method  used  hy  Benham,  Dowden  and  Kunzman  to  predict  the 
critical  pressure  and  temperature  of  a  naturally  occurring  hydrocarbon  mixture 
is  adequate. 

3)  If  a  more  precise  calculation  of  critical  temperature  and 
pressure  is  required,  either  Organick's  or  Eilerts'  procedure  should  be  used. 

4)  The  most  important  single  factor  to  be  considered  when 
calculating  the  critical  conditions  in  a  porous  media  is  the  composition  of 
the  fluid  at  any  particular  time. 

5)  The  conf  irmation  of  miscibility  conditions  by  laboratory  core 
testing  necessitates  the  elimination  or  suppressment  of  all  causes  of  poor  core 
sweep  efficiency.  The  observation  of  low  recoveries  is  not,  by  itself,  con- 
elusive  evidence  of  lack  of  miscibility. 

6)  At  viscosity  ratios  above  10.0,  the  mechanism  of  displacement 
in  a  sandstone  system  is  primarily  controlled  by  the  presence  of  viscous  fingers. 
Thus  any  attempt  to  correlate  a  "length  of  mixing  zone"  is  meaningless. 
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7)  At  viscosity  ratios  less  than  10.0  the  mechanism  of 
displacement  more  closely  resembles  the  ideal  piston  like  situation. 

8)  The  oil  recovery  to  breakthrough  of  the  displacing  phase 
at  high  viscosity  ratios  is  similar  in  magnitude  to  that  reported  for  water 
floods  with  the  same  viscosity  ratio.  However  the  subordinate  production  of  a 
miscible  flood  is  more  efficient  than  that  of  a  water  flood. 

9)  From  a  practical  point  of  view,  any  efforts  placed  on  the 
elimination  of  gravitational  effects  and  viscous  fingering  in  field  tests 
should  result  in  more  favorable  recoveries  by  miscible  displacement. 


RECOMMEND  AT  IONS 


Because  of  the  increasing  importance  of  recovering  more  oil  from 
known  reservoirs,  a  continuing  program  investigating  the  limitations  and 
advantages  of  a  miscible  type  system  is  justified.  Although  the  economics 
of  a  secondary  recovery  project  involving  miscible  displacement  are  sometimes 
questionable,  the  efficiency  of  an  ideal  miscible  flood  is  certainly 
attractive.  With  this  in  mind  it  is  recommended  that  the  research  program 
concerned  with  miscible  displacement  be  continued  and  broadened  to  include 
some  of  the  following  aspects:- 

1)  an  intensified  investigation  of  the  mechanism  of  displacement 
in  the  region  of  viscosity  ratios  less  than  10.0  using  a  vertical  reservoir 
and  additional  pressure  transducers. 

2)  an  investigation  of  the  results  achieved  using  enriched  gas 
mixtures  as  the  displacing  phase. 

3)  a  formulation  of  a  correlation  relating  pressure  drop  in  the 
porous  media  to  composition  of  fluid  flowing  in  the  system. 

4)  additional  investigation  of  the  growth  of  mixing  zones  at 
various  displacement  rates  employing  longer  cores. 
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NOMENCLATURE 


B 

n' 

Bi 

W 

ra. 


M- 


Ei 


'  cm 


T 

cm 

Aij 

Tc 

To. 


n 


pc 


n 


mT. 


mP 


c 
M 

M* 

T 

xce 

Pce 

% 

Bp 

N 


molal  average  boiling  point,  R 
mole  fraction  of  component  "i" 

o 

normal  boiling  point  for  component  " i" ,  R 

weight  average  equivalent  molecular  weight 

weight  fraction  of  component  "i” 

equivalent  molecular  weight  of  component  "i" 

measured  critical  pressure,  psia 

o 

measured  critical  temperature,  R 

characterization  factor  for  binary  mixture  of  components  "i"  and  "j* 

o 

true  critical  temperature,  R 
true  critical  pressure,  psia 

o 

molal  average  critical  temperature,  R 
molal  average  critical  pressure,  psia 

o 

weight  average  critical  temperature,  R 
weight  average  critical  pressure,  psia 
molecular  weight 
arbitrary  variable 

o 

estimated  critical  temperature,  R 
estimated  critical  pressxire,  psia 

normal  boiling  point  of  heavy  ( liquid ) portion  of  mixture,  R 

normal  boiling  point  of  light  (vapor)  portion  of  mixture,  °R 

o 

normal  boiling  point  of  entire  mixture,  R 

.  o 

mole  fraction  of  mixture  in  vapor  phase  at  14.7  psia  and  75  F 


T 


cc 


calculated  critical  temperature,  R 
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P=c 

= 

calculated  critical  pressure,  psia 

te 

SB 

0 

equivalent  critical,  temperature,  R 

A. 

S3 

mole  fraction  of  component  in  liquid  phase 

Q 

IS 

injection  rate,  cc/hr 

^0 

cs 

viscosity  of  oil,  cps 

^S 

S3  ' 

viscosity  of  solvent,  cps 

k 

» 

permeability,  md 

P 

up 

a 

upstream  core  pressure,  psig 

pdown 

= 

downstream  core  pressure,  psig 

L 

= 

length  of  core,  cm. 

A 

s 

2 

cross-sectional  area  of  core,  cm 

GOR 

a 

gas -oil  ratio,  cc/cc 

see 

= 

standard  cubic  centimeters  at  1^.7  psig  and  60  F 

stcc 

9 

stock  tank  cubic  centimeters  at  atmospheric  conditions 
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CALCULATION  OF  CRITICALS  BY  ORGANICK'S  METHOD  FOR  C- 
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TABLE  4 -a 


CALCULATION  OF  CRITXCALS  BY  MODIFIED  KURATA-KATZ  METHOD  FOR  C-l 


(1) 

(2) 

(3) 

w 

(5) 

(6) 

(7) 

(8) 

Component 

Mol.  Fr. 

Mol.  Wt. 

Tc 

(°R) 

Pc 

(psia) 

Mol.  Av. 
Mol.  Wt. 

=  (2)  x  (3) 

Ten 

=  (2)  x  (4) 
°R 

Pen 

=  (2)  X  5 

psia 

C°2 

0.0035 

44.01 

547.5 

1070.2 

0.154 

1.91 

3.74 

Cl 

0.6433 

16. 04 

343.2 

673.1 

10.300 

220.00 

425.00 

c2 

0.0638 

30.1 

549.9 

708.3 

1.920 

35.00 

45.20 

c3 

0.0605 

44.1 

666.0 

617.4 

2.670 

40.40 

37.40 

iC4 

0.0107 

58.1 

734.6 

529.1 

0.621 

7.85 

5.65 

nCL 

0.0458 

58.1 

765.7 

550.1 

2.660 

35.10 

25.20 

iCc 

O.OI89 

72.2 

829.6 

483.5 

1.360 

15.65 

9.14 

nC5 

0.0215 

72.2 

.846.2 

489.8 

1.550 

18.21 

10.50 

c6 

0.0366 

86.2 

914.2 

440.1 

3.160 

33.40 

16.10 

c7+ 

0.0954 

ll4.0^a) 

1010^) 

380. o( 

c)l0.86o 

96.20 

36.20 

(a)  Obtained  from  Table  1- 

(b)  and  (c)  from  Figure  30. 

,  ngaa^5) 

using  S. 

35.255 

G.  =  0.71 

503.72 

614.13 

(l)  Te  from  Fig.  5?A,  NGAA  =  595  R 


(2)  Tfi  .  595  .  1<l82 


Ten  503.72 


(3)  Tf  (r) 

m—  corrected  from  Fig.  62,  NGAAV  ' 
xcn 

f  c;  %  O 

(4)  True  Critical  Temp,  from  Fig.  6l  NGAA^ " '  =  6l6  R 


(5)  From  Fig.  4,  Reference 


(8) 


-ZsJLZsn-  =  i4o 


(6)  P, 


cc 


616  X  614.13 

140 


2700  p.s.i,a. 
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TABLE . 5 -A 

CALCULATION  OF  CRIEICALS  BY  MODIFIED  KORA TA-KATZ  METHOD  FOR  CH~P°843 


Component 

Mol.  Fr. 

Mol.  Wt. 

Tc 

(°E) 

Pc 

(psia) 

Mol.  Av. 
Mol.  Wt. 

-  (2)TxD(4) 

p 

-  (2)  x  (5) 

COp 

0.00794 

44.01 

547.5 

1070.2 

0.350 

4.347 

8.497 

n2 

0.01375 

28.02 

227.0 

492.2 

0.387 

3.121 

6.768 

Cl 

0.76432 

16. 04 

343.2 

673.1 

12.250 

262.315 

514.464 

C2 

0.07923 

30.1 

549.9 

708.3 

2.385 

43.569 

56.119 

c3 

0.04301 

44.1 

666.0 

617.4 

1.900 

26.645 

26.554 

iC4 

0.01198 

58.1 

734.6 

529.1 

0.706 

8.801 

6.339 

nC4 

0.01862 

58.1 

765.7 

550.1 

1.100 

14.257 

10.243 

iC5 

0.00937 

72.2 

829.6 

483.5 

0. 676 

7.773 

4.530 

nCr 

O.OO78I 

72.2 

846.2 

489.8 

0.564 

6.609 

3.825 

c6 

0.01405 

86.2, 

914.2 

440.1 

1.210 

12.845 

6.183 

c7  + 

0.02992 

127. 3V a) 

106o(t>) 

360(c) 

3.820 

31.715 

10.771 

1.00000 

25.348 

421.997 

654.293 

(a)  Obtained  from  Table  1-A 


(b)  and  (c)  from  Figure  30»  N GAAV^  using  S.G.  *=  0.74 


(5) 


(1)  Te  from  Figure  57A,  NGAA^-’  =  492°R 

(2)  S&  =  1.167 


m 

cn 


corrected  from  Figure  62,  NGAA^'*  =  1.068 

^cn 

(4)  True  Critical  Temperature  Figure  6l,  NGAA^'  =  520°R 

(5)  From  Figure  4,  Reference^) 

JkjLl^GL.  =  160 


(6)  P, 


cc 


520  x  654 

160 


2130  p.s.i.a. 
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Liquid  Phase  Gas  Phase  Point 
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APPENDIX  B 
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TABLE  1-B 


DESCRIPTION  OF  BEREA  SANDSTONE  CORE 


Core 

Length 

= 

3 6  in.  =  $l.h  cm. 

Core 

Diameter 

= 

2  in.  -  5.18  cm. 

Core 

Cross-sectional  Area 

= 

21.05  em.^ 

Core 

Bulk  Volume 

= 

113  cu.  in.  =  1925.0  cc. 

Core 

Porosity  (Average) 

= 

16.30 

Core 

Permeability  to  Air 

= 

29 3-9  millidareies 

Core 

Permeability  to  Liquid  (varsol) 

= 

266  millidareies 
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TAHLE  2-B 


FLUID  PROPERTIES 


(1)  (2)  (3)  (4)  (5)  (6) 


Property 

Crade  Oil 
Sample 

Varsoi 

Pentane 

Propane 

Esso  P~40 

M.W. 

197*78  (calc) 

- 

72.146 

44.094 

- 

Density  at 

60/60 

0.848  (calc) 

- 

0.63116 

0. 5077 

0.8165 

(gm/cc) 

s> 

Density  at  76  F 

O.8536 

0.7809 

a 

0.8090 

( gm/cc ) 

Kinematic  p 

1.2200 

9.097 

«=> 

«=» 

1.969 

( cent i stokes) 

Absolute  p 

0.9528 

7.765 

0.23 

0.110  at  1500 

psia 

1.53^ 

(centipoise) 

Mol.  °lo  N 

0.0 

0.108  at  1000 

psia 

Mol.  fo  C02 

0.0 

- 

- 

- 

- 

Mol.  io  H2S 

0.0 

= 

- 

- 

■= 

Mol.  io  C± 

0.0 

~ 

- 

- 

Mol.  i  C2 

0.0 

- 

- 

- 

- 

Mol.  io  Co 

0.0 

- 

- 

Mol.  io  ICij. 

0.0 

- 

= 

- 

Mol.  io  nC^ 

6.60 

- 

- 

- 

Mol.  io  1C5 

3.W 

- 

- 

Mol.  io  nCc 

6.01 

- 

to 

- 

Mol.  io  Cg 

7.80 

- 

- 

- 

- 

Mol.  io  C7+ 

76.11 

- 

- 

- 

~ 

S.G.  C7+ 

0.877 

- 

- 

- 

M.W.  c|+ 

237.0 

“ 

“ 

- 
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AIR  PERMEABILITY  DATA 


IQ 

ft 

O 


UN 

CM 

CO 

CO 

H 

o 

o 


Sh 

"H 


CO 


U-\ 

o 


o 


CD 

u 

o 

o 

«H 

o 

ed 

<D 


ft 

o 


p 

II 

o 

rH 

K 

t- 

co* 


PQ 

i 

CO 


p 

rH 

on 

II 


o 

O 

«H 

o 


a 


t— 

CO 

X 

P 

o 

H 

on 

II 

P 


w 

CD 

“H 

O 


f-1 


•H 


H 

ri 


p  d  <d 

UN  Ijjj  ^ 

lop 

Qffi  II! 


O  O 

o  CD 
'  to 


m  to 

P  q  to 

CS3  0  0)  O  G 

_  P  H  o  CD 

q?  <:  i ±n  i  —  «j 


to  to 
P  O  to 
0}  0) 

'Har  is  & 


O.J 


CO 

p 

in 


to 

*3 


$ 

oJ 

o 

to 

•3 


tp 

rH 

tP 

co 

§ 

r“l 

un 

to 

P 

on 

p 

o\ 

-4 

-4 

un 

rH 

P 

<8 

CO 

OJ 

CO 

CO 

< — 1 

8 

e 

CM 

O 

o 

O 

d 

o 

d 

d 

d 

o 

d 

o 

O 

ITS 

o 

ir>. 

o 

to. 

O 

o 

UN 

© 

9 

© 

• 

• 

© 

• 

« 

o 

O 

o 

o\ 

oo 

VO 

rH 

ON 

o 

CO 

4 

on 

CO 

on 

O', 

on 

On 

00 

CO 

o 

O 

CM 

Oi 

CVS 

OJ 

CVS 

CM 

CM 

C-J 

ro 

CO 

5$ 

IP 

CO 

4 

vo 

CO 

vo 

LfN 

t*— 

CO 

CM 

ON 

4 

Ed 

4 

CO 

ON 

vo 

vo 

CO 

CO 

LT\ 

o\ 

4 

vo 

O 

p 

vo 

• 

o 

O 

• 

• 

© 

• 

• 

o 

• 

o 

H 

rH 

rH 

CM 

CM 

CM 

CO 

co 

co 

4 

LfN 

O', 

vo 

O 

4 

LfN 

rH 

LfN 

p 

C\ 

<s 

rH 

CM 

ro 

rH 

CO 

rH 

CO 

o 

oo 

CO 

O 

CM 

-4* 

vo 

O 

p 

C'— 

o 

« 

• 

© 

© 

• 

© 

• 

• 

• 

1 — 1 

H 

« — 1 

CM 

CM 

0J 

CM 

CO 

CO 

CO 

CM 

LfN 

o 

UN 

co 

O 

-4 

O 

o 

LfN 

IP 

o 

co 

O 

P 

On 

P 

LfN 

CO 

© 

• 

• 

» 

© 

• 

• 

• 

• 

o 

4 

CO 

o 

to 

vo 

CO 

ON 

0J 

vo 

kr 

H 

H 

H 

H 

rH 

CJ 

CM 

ON 

UN 

ON 

o 

UN 

O 

00 

vo 

co 

On 

co 

H 

CO 

CO 

• 

© 

o 

• 

• 

» 

0 

o 

• 

• 

-4 

00 

rH 

LfN 

ON 

CO 

co 

UN 

tP 

vo 

H 

H 

H 

CM 

CM 

CO 

p 

UN 

O 


UN 

CO 

ON 

o 

5 

& 

ON 

LfN 

o 

IP 

o 

rH 

UN 

ON 

O 

VO 

o 

rH 

<o 

4 

UN 

VO 

00* 

ON 

CO 

VO 

d 

H 

rH 

CM 

CM 

cB 

8 

P 

OJ 

O 

CM 

CO 

H 

o 

UN 

UN 

CJ 

ON 

8 

VO 

o 

t— 

ON 

ON 

o 

O 

H 

rH 

P 

co 

• 

o 

© 

• 

• 

© 

© 

• 

• 

• 

o 

O 

rH 

1 - 1 

H 

rH 

rH 

1 — 1 

rH 

rH 

vo 

CO 

00 

00 

CO 

co 

co 

oo 

oo 

00 

CO 

H 

rH 

ON 

CO 

CO 

t- 

UN 

VO 

CO 

00 

CO 

vo 

ON 

co 

vo 

On 

UN 

H 

UN 

• 

• 

• 

• 

• 

• 

© 

o 

© 

• 

H 

CM 

CM 

CM 

co 

CO 

co 

-4 

UN 

UN 

Page  B-3 


Average  Permeability  to  Air  =  293.9  md. 


4-b  liquid  (varsol)  permeability  data 
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Average  Permeability  to  Liquid  (Varsol)  =  2 66„4  md. 


TABLE  5-B  PRESSURE  PROFILE  DATA  FLOWING  VARSOL 
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PRESSURE  PROFILE  DATA  FLOWING  VAR  SOL 
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TABLE  7=B  PRESSURE.  PROFILE  DATA  FLOWING  PROPANE 
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Length  of  Core  =  36  inches 
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TABLE  9-B  PRODUCTION  DATA 
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13-B  PRODUCTION  DATA 
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Recovery  of  oil  after  displacement  of  1,5  pore  volumes  =  "(Ofo  0IP 


TABLE  l4-B 


PRESSURE  PROFILE  DATA 


Displacement  No.  D 


Propane  Displacing  Crude  Oil  Q 


60  cc/hr 


Length  of  Core  =  36  in. 

Sec.  1=6  in.  Sec.  2  =  12  In.  Sec.  3  =  12  in. 


Sec.  4=6  in. 


(1)  (2)  (3)  (b)  (5)  (6) 


Time 

Pressure  Differential. 

per  inch  (AP/AL) 

Sec. 

1  Sec.  2 

Sec.  3 

Sec.  4 

Total 

(min; sec) 

(  P. 

s.i.d./in. ) 

0:00 

1.13 

1.53 

0.96 

1.23 

1.22 

4;  30 

O.83 

1.52 

1.00 

1.20 

1.20 

7:00 

0.60 

1.52 

1.00 

1.20 

1.14 

10:00 

0.40 

1.19 

0.99 

1.18 

0.98 

13:45 

0.30 

1.05 

0.99 

1.18 

0.93 

15:45 

0.30 

0.99 

0.99 

1.18 

0.91 

17:45 

0.30 

0.92 

0.98 

1.18 

0.88 

19:45 

0.30 

'  O.83 

O.98 

1.17 

O.85 

22:00 

0.28 

0.77 

O.98 

1.17 

O.83 

25:20 

O.30 

0.68 

O.98 

1.17 

0.80 

29:45 

0.33 

0.53 

O.98 

1.17 

0.75 

32:25 

O.38 

0.45 

O.98 

1.17 

0.74 

37:00 

0.52 

0.30 

0.93 

1.17 

0.69 

42:30 

0.60 

0.22 

0.82 

1.17 

0.64 

49.45 

0.57 

0.17 

O.67 

1.17 

0.57 

56:15 

0.57 

0.l4 

0.52 

1.15 

0.51 

60:00 

0.53 

0.13 

0.43 

1.13 

0.47 

64: 40 

0.50 

0.12 

0.37 

1.08 

0.43 

68:35 

0.48 

0.12 

0.32 

0.98 

0.39 

71:20 

0.47 

0.11 

0.28 

0.90 

0.36 

75:10 

0.47 

0.10 

0.26 

0.77 

0.33 

79:15 

0.45 

0.09 

0.23 

O.58 

0.28 

81:50 

0.45 

0.09 

0.22 

0.52 

0.26 

84:30 

0.45 

0.09 

0.21 

0.47 

0.25 

85:30 

0.45 

0.08 

0.20 

0.45 

0.24 

91:15 

0.42 

0.08 

0.17 

O.38 

0.21 

96:15 

0.42 

0.08 

0.16 

0.35 

0.20 

99:35 

o.4o 

0.08 

0.16 

0.32 

0.20 

100:15 

0.40 

0.08 

0.16 

0.32 

0.20 

102:30 

0.40 

0.07 

0.15 

0.30 

0.19 

105:30 

o.4o 

0.07 

0.l4 

0.28 

0.18 

108: 30 

0.40 

0.07 

0.14 

0.28 

0.18 

113:20 

0.40 

0.07 

0.13 

0.25 

0.17 

116:10 

0.40 

0.07 

0.13 

0.25 

0.17 

119: 25 

0.38 

0.06 

0.12 

0.23 

0.16 

122:40 

0.38 

0.06 

0.11 

0.22 

0.16 

125:45 

0.38 

0.0  6 

0.11 

0.22 

0.16 
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TABLE  14-B  ( Continued ) 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

Time 

Pressure  Differential  per  in 

(AP/AL) 

Sec.  1 

Sec.  2 

Sec.  3 

Sec.  4 

Total 

(min: sec) . 

'  '  l 

( 

p.s.i.d./in. 

) 

130: 30 

0.37 

0.05 

0.10 

0.20 

o.i4 

134:40 

0.37 

0.05 

0.10 

0.20 

0.l4 

137-02 

0.37 

0.05 

0.10 

0.20 

0.14 

l4o : 19 

0.37 

0.05 

0.09 

0.18 

0.14 

145:29 

0.37 

0.05 

0.08 

0.18 

0.14 

149:07 

0.35 

0.05 

0.08 

0.18 

0.13 

153:09 

0.35 

0.05 

0.08 

0.18 

0.13 

156:45 

0.35 

0.05 

0.08 

0.18 

0.13 

160 : 18 

0.33 
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0.13 

163:43 

0.33 

0.05 

0.07 

0.17 

0.12 

169:07 

0.33 

0.05 

0.07 

0.17 

0.12 

173:03 

0.33 

0.05 

0.07 

0.17 

0.12 

180:30 

0.33 

0.04 

0.06 

0.15 

0.11 

184:45 

0.33 

0.04 

0.06 

0.15 

0.11 
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0.32 
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0.06 

0.12 

0.11 
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0.32 

0.04 

0.06 

0.12 

0.11 

208:10 

0.30 

0.04 

0.05 

0.10 

0.10 

222:20 

0.28 

0.04 

0.04 

0.08 

0.09 

236:00 

0.28 

0.04 

0.04 

0.08 

0.09 

251:00 

0.28 

0.04 

o.o4 

0.08 

0.09 

Pressures  Became  Erratic 
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15-B  PRESSURE  PROFILE  DATA 
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16-B  pressure  profile  data 
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